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1.  Overview  of  Activities  and  Progress 
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to  expand  the  scope  of  activity.  Additional  research  and  development  tasks  were  undertaken,  as 
follows. 


C Provision  of  Revised  BATS  Parameters  from  Field  Data 


(b)  Shortcomings  in  the  representation  of  heterogeneous  vegetation  cover  using  aggregation  rules 
are  possible  for  landscapes  and  climates  in  which  there  are  either  large  differences  in  soil  moisture 
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(ii)  Aggregation  rules  are  relatively  straightforward  in  the  case  of  patch-scale  heterogeneity  of 
vegetation-dependent  controls  on  surface  exchanges,  but  aggregation  of  soil  hydraulic 
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3.  List  of  Publications  Resulting  from  this  Contract 
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University  of  Arizona,  Tucson,  150p. 
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(Copies  of  papers  2-17  are  provided  herewith). 
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crop  were  selected  for  the  aggregation  studies.  The  parameter  values  associated  with 
the  selected  land  covers  are  given  in  Table  1. 

Soil  texture  class  7 and  soil  color  class  5 were  selected  as  being  appropriate  for  the 
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Table  3 

BATS  state  variables,  driving  variables,  and  parameters 
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land  surface  models.  The  the  surface  water  balance  of  these  cover  classes  was  calculated  conservatively); 

t simple  aggregation  rules,  acceptable  behavior  was  always  observed  for  mixes  of  short  and  tall  grass  and 

>1,  were  successful  at  pro-  agricultural  crops.  We  suspect  that  this  is  likely  a more  general  result,  especially 

variety  of  realistic  hydro-  when  seeking  aggregate  descriptions  of  the  restricted  mixture  of  land  covers  which* 
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and  surface  models.  The  the  surface  water  balance  of  these  cover  classes  was  calculated  conservatively); 

simple  aggregation  rules,  acceptable  behavior  was  always  observed  for  mixes  of  short  and  tall  grass  and 

*1,  were  successful  at  pro-  agricultural  crops.  We  suspect  that  this  is  likely  a more  general  result,  especially 

variety  of  realistic  hydro-  when  seeking  aggregate  descriptions  of  the  restricted  mixture  of  land  covers  which 
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MAPPING  SURFACE  COVER  PARAMETERS  USING  AGGREGATION 


RULES  AND  REMOTELY  SENSED  COVER  CLASSES 


Abstract 


CM 


1.  Introduction 


» 


schemes  actually  in  use  within  General  Circulation  Models  (GCMs). 


The  Biosphere- Atmosphere  Transfer  Scheme  (BATS,  Dickinson  et  al.,  1986)  is  the  land-surface 
parameterization  used  in  the  National  Center  for  Atmospheric  Research’s  (NCAR’s)  Community 
Climate  Model.  In  normal  application  of  BATS,  the  values  of  vegetation  parameters  are  based  on 
the  single  dominant  surface-cover  in  each  grid  cell.  A set  of  hypothetical  aggregation  rules  for 


a variety  of  hydrometeorological  conditions  at  this  site.  However,  they  did  not  work  well  when 
artificially  wet  (irrigated)  patches  were  combined  with  naturally  dry  soil  patches. 


In  the  Arain  et  al.  (1996)  study,  the  application  of  aggregation  rules  for  some  important  BATS 
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On  the  basis  of  the  success  reported  in  the  current  (and  previous)  study,  the  BATS-specific 


aggregation  rules  were  then  applied  to  the  1 km  land-cover  data  set  of  Loveland  et  al.  (1991)  to 
calculate  grid-average  values  of  BATS  vegetation  parameters  for  2.8  x 2.8  and  lxl  grid 
squares  across  the  conterminous  United  States  (U.S).  A comparison  of  the  ensuing  aggregate 


times  when  some  of  the  meteorological  and  soil  moisture  data  were  missing  at  one  or  both  sites 


2.2  Validation  of  BATS  at  the  Amazon  Sites 


necessary  to  initiate  nine  state  variables  (see  Table  4).  Surface  soil  temperature,  deep  soil 


5a).  In  this  study,  the  philosophy  adopted  was  to  make  the  minimum  parameter  changes 
necessary  to  achieve  acceptable  agreement  between  observations  and  model  calculations.  In 
practice,  the  simulations  of  the  pasture  energy  fluxes  were  greatly  improved  merely  by  changing 


moisture  and  meteorological  forcing  data  were  available  at  the  Reserva  Ducke  forest  site. 
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The  interface  height  between  BATS  and  the  atmospheric  boundary-layer  model  was  set  to  45  m. 
[For  a detailed  discussion  of  the  sensitivity  of  simulated  fluxes  to  the  choice  of  this  model 


interface  height,  see  Arain  (1994)  and  Arain  et  al.  (1996)].  Initial  vertical  profiles  of  potential 


site,  but  rainfall  is  seasonal.  By  chance,  one  dry  period  during  the  1990  data 
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nding  aggregate  vegetation  cover  also  did  not  show  evidence 


point,  while  there  was  still  plenty  of  water  available  for  evaporation  in  the  deeper-rooted  forest 
patch.  The  variation  in  sensible  heat,  latent  heat,  soil  heat,  and  momentum  flux  across  the 
modeled  domain  from  two-patch  and  homogeneous  cover  is  shown  in  Figure  2.  The  differences 


vegetation  parameters  were  used  (see  Section  2.2).  This  figure  shows  that  the  aggregation  rules 


orked  quite  well  when  neither  of  the  vegetation  patches  has  soil  moisture  stress,  but  agreement 


In  order  to  investigate  the  sensitivity  of  aggregation  rules  to  large  differences  in  soil  moisture  in 
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stomatal  resistance  for  short  grass  (40  sm’1)  and  evergreen  forest  (150  sm  ) gives  the  values  of 


days  and  three  times  of  the  day  for  which  tests  were  made.  On  the  basis  of  this  evidence,  a linear 
aggregation  rule  is  recommended  for  this  BATS  parameter,  and  linear  aggregation  of  minimum 
stomatal  resistance  was  adopted  and  applied  for  the  remainder  of  the  analysis  described  below. 


3.  Application  of  Aggregation  Rules  for  USA 


landscape,  and  they  arise  because  the  root  resistance  assumed  in  BATS  is  a non-linear  function  of 
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questionable  in  two  specific  cases,  namely  for  landscapes  and  climates  in  which  there  are  either 
(a)  large  differences  in  soil  moisture  and  with  wet  and  dry  soils  covering  approximately  equal 
areas  and  with  soil  moisture  differences  greater  than  10-20%  of  the  rooting  zone  soil  moisture 
store;  or  (b)  approximately  equal  areas  of  shallow-rooted  and  deep-rooted  vegetation  growing  in 
a climate  where  the  shallower-rooted  vegetation  is  under  soil  moisture  stress. 


We  obtained  aggregate  BATS  parameter  values  for  2.8°  GCM  grids  by  computing  the  area- 
average  of  the  1 km  BATS  values  obtained  from  the  USGS-EROS  land  cover  data  base  described 
above.  Aggregate  values  of  vegetation  parameters  (Table  3)  were  computed  using  the 


aggregation  rules  given  in  Table  1 (but  assuming  a linear  aggregation  rule  for  minimum  stomatal 


total  number  of  (1  km  x 1 km)  pixels  for  each  cover  class.  In  the  case  of  cover  types  assigned  a 
number  20-28  in  the  USGS-EROS  data,  50%  of  the  area  of  each  pixel  was  assigned  to  be  one 
cover  type  and  50%  to  the  other  cover  type  in  these  calculations.  The  number  of  remotely  sensed 


aggregate-cover  using  USGS-EROS  data  is  illustrated  in  Figure 


vegetation  cover  is  influenced  little  by  using  these  different  approaches  of  parameterization  for 


The  calculated  values  of  area-average  albedo  and  warm  season  leaf  area  index  (Figures  8c-d)  also 
show  marked  differences,  especially  in  the  western  United  States.  The  use  of  default  cover  class 
assigned  in  CCM2,  which  is  often  short  grass  or  shrub  land,  neglects  the  presence  of  other 


surface  cover  types  such  as  desert  and  evergreen  forest.  However,  the  presence  of  these  covers  is 
reflected  in  the  parameter  values  calculated  for  aggregate  cover.  The  assigned  values  of  cold 
season  leaf  area  index  and  aerodynamic  roughness  length  exhibit  marked  differences  (of  up  to 


between  parameters  in  the  case  of  this  smaller  1°  x 1°  grid  mesh.  The  most  noticeable  distinction 
between  Figure  8 and  Figure  9 occurs  in  the  western  United  States  for  vegetation  fraction,  leaf 
area  index,  albedo,  and  minimum  surface  resistance.  This  distinction  is  not  related  to  the  use  of 
aggregate  covers,  but  is  due  to  a mismatch  between  the  single  cover  currently  assigned  within 
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In  BATS,  air  temperature,  dew  point  temperature,  surface  pressure,  wind  velocity,  precipitation, 
and  downward  short-wave  and  long-wave  radiation  are  used  as  forcing  variables,  and  these 
variables  are  available  as  6-hourly  time  series  on  the  ISLSCP  CDROM.  They  were  interpolated 


to  provide  hourly-mean  values  using  a cubic  spline  interpolation  method.  The  interpolated 


parameters  corresponding  to  these  soil  classes  was  then  made  following  Dickinson  et  al.  (1986). 
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heat  fluxes  (Figures  lOg-h)  are  calculated  to  be  negative  in  the  eastern  USA  but  are  positive  and 
approximately  equal  to  the  latent  heat  flux  in  the  western  USA.  The  BATS  simulated  fluxes  were 
also  compared  with  estimated  surface  fluxes  provided  on  the  ISLSCP  CDROM  which  were 


extracted  from  the  European  Center  for  Medium-Range  Weather  Forecast  (ECMWF)  operational 


the  page  (Figures  1 la-b)  are  for  latent  heat,  and  those  on  the  right  side  of  the  page  (Figures  1 lc- 
d)  are  for  sensible  heat.  The  difference  in  calculated  fluxes  is  within  ±10%  in  many  locations  and 
in  both  seasons.  However  there  are  significant  differences  in  the  sensible  heat  flux  in  the  summer 


season  (Figure  11c),  which  are  particularly  noticeable  in  parts  of  the  western  USA  and  are 
consistent  with  the  large  difference  in  surface  albedo  and  leaf  area  index  for  the  two 
parameterizations  in  this  region.  There  is  also  an  appreciable  shift  in  the  Bowen  ratio  in  some 


little  difference  throughout  the  entire  two-year  period,  but  in  Arizona  including  the  influence  of 


minimum  stomata!  resistance.  Aerodynamic  roughness  length  is  the  exception:  in  this  case. 


logarithmic  aggregation  rule  is  required  (Table  1).  Experiments  conducted  for  the  ABRACOS 
site  were  effective  tests  of  aggregation  rules  for  stomatal  resistance,  aerodynamic  roughness 


Tfr 


<N 


but  the  proportion  of  irrigated  land  within  a GCM  grid  is  usually  small  and  its  influence  on  the 
area-average  fluxes  likewise.  Large  spatial  variability  in  plant-available  moisture  might  have 
natural  origin,  perhaps  resulting  from  intermittent,  spatially  heterogeneous  rainfall.  A parallel 
study  (White  et  al.,  1996)  directly  addressed  this  issue  for  the  extreme  case  of  convective  rain 
falling  in  a semi-arid  environment.  White’s  study  demonstrated  that  sufficiently  large  differences 


extreme  situations  it  is  perhaps  necessary  to  find  a way  to  aggregate  soil  moisture  and  vegetation 


Further  research  is  required  to  evaluate  whether  the  difference  in  surface  fluxes  given  by  using 
aggregate  (as  opposed  to  most-common  cover)  parameters  has  a significant  influence  on  modeled 
climates  given  by  GCMs.  However,  global  application  of  aggregation  rules  to  improve  GCM 
performance  (in  conjunction  with  upcoming  remotely  sensed  AVHRR  Path-Finder  and  EOS 


MODIS  data  sets)  is  at  this  stage  arguably  more  profoundly  limited  by  shortcomings  in  the 


acquired  and  distributed  the  USGS-EROS  land  cover  and  the  ABRACOS  data  sets 
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Figure  6.  Comparison  of  the  calculated  fluxes  given  with  three  alternative 
aggregation  rules  for  minimum  stomatal  resistance  on  (a)  October  30, 
1990  (a  day  with  very  dry  soils),  and  (b)  August  28,  1991  (a  day  with  wet 
soils).  Calculated  latent  heats  are  shown  with  the  lines  without  symbols, 
and  sensible  heat  fluxes  are  shown  with  lines  with  crosses. 
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Figure  12  Monthly  average  latent  heat  and  sensible  heat  fluxes  for  1987  and  1988 
calculated  with  a stand-alone  BATS  model  forced  with  ISLSCP  CDROM 
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Table  6.  Days  and  times  for  which  aggregate  and  patch  simulations  were  performed  and 
observed  soil  moisture  at  corresponding  mid-day  time. 


Executive  Summary  of  the  Tucson  Aggregation  Workshop 


efficient  use  of  highly  resolved  spatial  data,  or  by  the  desire  to  proceed  confidently  without 
utilization  of  detailed  data.  In  other  words,  it  seeks  to  address  the  question,  "How  can  we  model 
spatially  variable  processes  using  a grid  size  which  is  coarse  enough  to  be  economical,  yet  fine 


h that  results  are  not  affected  by  sub-grid  scale  variability?"  However,  the  topic  of 


and  the  International  Satellite  Land  Surface  Climatology  Project  (ISLSCP)  of  the  World  Climate 
Research  Program.  In  some  part,  the  organization  of  the  workshop,  and  in  major  part  the 
subsequent  publication  of  this  Special  Issue  was  also  supported  by  the  U.S.  National  Aeronautic 


to  these.  In  other  cases,  the  authors  chose  to  respond  to  the  audience's  questions  by  addressing 
these  in  the  main  body  of  the  paper  they  subsequently  prepared.  It  is  our  opinion  that  the  ensuing 
suite  of  papers  comprise  a focused  but  broad  overview  of  our  current  understanding  of  issues 
related  to  the  aggregate  representation  of  land-atmosphere  interactions,  and  it  is  with  satisfaction 
and  some  pride  that  we  propagate  them  in  this  Special  Issue  of  the  Journal  of  Hydrology. 
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The  purpose  of  many  aggregation  studies  is  to  provide  information  which  is  intended  to 


Remotely  sensed  vegetation  indices  contain  useful  information  on  the  bulk  stomatal 
resistance  and  photosynthetic  uptake  of  vegetation  (Sellers  et  al.,  1992),  but  the  roles  of 
vegetation  type  and  nutrition  on  the  interpretation  of  these  indices  requires  further 


Aggregation  of  remotely  sensed  measurements  in  sparse  canopies  can  be  accomplished 
ti  little  error  in  some  circumstances  (such  as  aggregation  of  surface  temperature  from  1 m2 


Summary  of  Progress  in  Aggregation  Research  and  Recommendations  for  Future 


There  has  been  significant  and  satisfactory  progress  in  developing  aggregation  techniques 


to  deal  with  vegetation-dependent  surface  properties  such  as  albedo,  aerodynamic  resistance, 


00 


sparse  and  dense  canopies  for  variables  derived  from  remotely  sensed  data.  W e do,  however. 
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Most  participants  favored  a "broad-front"  approach  combining  bold  application  of  current 


aggregation  strategies  at  global  and  regional  scales  and  continuing  mechanistic  research  into 
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the  soil  surfaces  in  the  model.  Both  models  describe  momentum  transfer  tradi 
in  terms  of  aerodynamic  roughness  which  is  prescribed  for  each  cover  type  in  B/ 
calculated  during  a pre-processing  phase  from  prescribed  canopy  morphology  in 
of  the  SiB  model. 


lentum  transfer  traditionally,  moisture  content  within  a rooting  layer.  Other  models,  such  as  that  of  the  French 

:ach  cover  type  in  BATS,  but  Meteorological  Service  (Noilhan  and  Planton,  1989),  are  similar,  but  parameterize  surface 

opy  morphology  in  the  case  resistance  in  a way  broadly  equivalent  to  BATS. 
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leling  studies  in  which  they  used  a I contamination  by  atmospheric  processes  for  example  - but  these  are  not  the  subject 
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Hydrological  Models,  Regional  Evaporation,  and  Remote  Sensing:  Let’s  Start  Simple 
and  Maintain  Perspective 
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describing  area-average  surface  exchanges  in  GCMs  is  arguably  less  than  that  required 
for  interpretative  hydrological  modeling,  because  land  surface  hydrology  is  just  one  of 
several  influences  on  climate,  and  GCMs  must  retain  an  appropriate  balance  in  the  way 
they  include  these  influences.  In  this  case,  an  appropriately  defined  averaging  procedure 
applied  to  component  land  covers  in  each  grid  square  might  suffice  to  define  the 
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one  hand,  irregular  topographically  determined  catchments  on  the  other)  can  be 
reconciled,  leaving  each  type  of  model  free  to  develop  in  its  own  optimum  and  balanced 


2.  The  Objectives  of  Hydrometeorological  Modeling 


° i 

e a 

0 V 

1 I 
i ! 

O O 

: ? 

6 

I 

E 


v 

8 

*d 


e 


l 


V 

5 


4) 

*d 
O 

6 

1 

a 

o 

u 

|>  ^ 
43  * 


% 


x "3 
^ * 


c 

o 

*^3 

§ «p4 

*d 

d 

H 


d 

§ 

•§ 

E 


1 


c 
o 

v>  (X 

5 

3 

-d 


Vi 

g 


*d 

i 

a 

^6 


d 


i 


3 

bi> 

C 

u 

•d 

o 

6 


n 


ro 


c 

o 

•d 

'■a 

> 

*-* 

o 

0) 


§ 


i 

a 

o 

a 

8* 

9 


u 

U 

a 


9 

i 

o 

u 

8 


*Eb 

o 

u 

£ 

00 

a 

■H* 


I 

a 

o 

.S 

u 
* ^ 

a 

jo 

rt 

V) 


d 

5 

V 

3 

8 

V 

.o 


U 

O 


w 

V 


So 

a 

u 

•a 

o 


I 

■d 

v 

ts 

• H 

"a 

E 

• *-* 
V) 

4> 

•d 


s 

0 

1 

c 

0 

1 


£ 

'$ 

13 

*d 

o 

*-» 

r 


I 

I 


•d 

4) 

N 


*Eh 

0 

Ml 

1 

I s 

5 p 

o 

a u 

.r  v 

I I 

5 .a 


1 „ 

i | 

Q.  « 


■8 

2 

* 


a 

•6  ° 

■a  3 

n 

8 J 

1 o 

•d 

d 


T3 

'8 

8* 

4) 

•s 

o 

JO 

c 

4) 

d 

* 

d 

JO 


o 

£ 

o 

u 

d 


to 


3 
8 
4) 

is 

-d 
4) 

b 

4) 

TJ 

! 

s f •* 

„ 1 1 

4)  o 

-gate 

§ T3  § 

I •&  i 

•a  1 U 

0 -3  > 

iii 

1 .i 

<«  U 

0 2 

1 1 

v 9* 

e § 

8 8 

| ? 

| I 

^ s 

•2  {j 


xi  -g 


r 


a 

o> 


c 

o 

'M 

Wl 


§ g 9 a 


8 


1 

J3 


« « 

J3  J3  Kh 

O >-  W> 

U U-l  A) 

a o w 


against  runoff  data,  the  required  hydrological  modeling  is  essentially  simple  and  large  scale 
and  that  there  is  no  real  need  to  couple  even  this  simple  routing  model  into  the  GCM. 
There  may,  however,  be  other  reasons  for  doing  this  coupling,  as  described  in  Sections 
23  and  2.4. 


A further  mode  of  validation  might  exist  in  which  an  advanced  hydrological  model  is 
used  to  synthesize  the  area  average  values  of  some  of  the  model-calculated  variables 
which  the  hydrological  model  shares  with  the  GCM.  Such  variables  might  include  soil 


Most  of  the  considerable  effort  which  has  so  far  been  directed  towards  improving  the 
representation  of  hydrology  in  GCMs  has  focused  on  providing  a more  realistic 
description  of  the  vertical  movement  of  water.  This  is  because  it  is  these  vertical  flows 
which  are  most  important  in  determining  the  net  return  and  the  rate  of  return  of  water 
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Similarly,  water  which  has  been  lost’  to  rivers  or  lakes  can  in  time  become  accessible 
to  the  atmosphere  again  if  the  level  of  the  water  body  drops  by  evaporation  to  reveal 
successively  deeper  layers  of  water.  This  process  could  be  important,  for  instance,  in  the 


case  of  many  shallow,  stagnant  lakes  in  high  latitude  North  America.  Otherwise,  stream 
water  drains  through  the  bed  of  the  river  or  lake  or  it  runs  off  down  river  out  of  the  grid 
square  and,  for  all  practical  purposes,  it  is  no  longer  available  to  support  evaporation 
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In  the  case  of  very  large,  flood-prone  rivers,  description  of  the 
intergrid  flow  mechanisms  which  allow  water  lost  to  river  flow  in 
one  grid  area  to  become  re-accessible  to  the  atmosphere  in  other 
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few  topography-determined  p; 


including  a few  ’representative  slopes',  perhaps  one  for  each  allowed  vegetation  cover, 
for  which  hillslope  flow  is  modeled  in  each  GCM  grid  square?  Progress  in  this  branch 
of  research  is  rapid,  and  it  is  likely  these  questions  will  be  answered  in  the  next  few 


significance  in  improving  the  prediction  of  global  environmental  change.  The  most 
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to  be  linear  functions  of  soil  water  content  and  because  a significant  proportion  of  those 
biogeochemical  releases  with  most  global  significance,  such  as  that  of  methane  from 
swamps  or  rice  paddies,  may  occur  from  areas  so  small  that  they  are,  in  global  terms, 
essentially  point  sources.  These  two  factors  argue  for  the  use  of  detailed  high-resolution 
hydrological  modeling  to  address  this  issue,  but  it  is  not  clear  that  implementing  models 


with  this  complexity  across  the  globe  is  justified,  especially  if  it  were  also  proposed  to 
include  such  a submodel  in  (GCM)  real-time.  Perhaps  some  compromise  will  emerge 
in  the  course  of  the  next  decade  involving  detailed  submodeling  for  key  source  areas,  but 
using  a model  time  step  which  is  less  than  that  used  in  the  atmospheric  component 
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Institute  (Shuttleworth,  1993). 


It  is  first  necessary  to  consider  which  aspects  of  the  energy-water  interactions  of 
vegetation-covered  continental  surfaces  are  truly  essential  to  include  in  SVAT  submodels 
for  the  purposes  of  adequate  climate  and  weather  prediction.  As  stated  earlier,  the  key 
characteristic  of  such  aspects  (and  the  parameters  which  determine  them)  is  that  they 
must  strongly  influence  the  extent  to  which  energy  and  water  is  available  to  the 


to  acknowledge  and  include  two  types  of  stores  in  models  of  the  land-atmospnere 
interaction.  The  first  is  the  energy  and  water  stored  in  comparatively  easily  accessible 
form  in,  or  in  the  case  of  water  on,  the  vegetation  and  soil  surfaces.  The  second  store 
is  the  energy  and  water  held  in  less  readily  accessible  form  in  the  soil,  but  only  to  the 
depth  where  they  are  still  (ultimately)  available  for  return  to  the  atmosphere. 


Because  the  timing  of  the  release  of  available  energy  and  water  is  important,  the  control 
exerted  by  the  vegetation  and  soil  in  determining  the  rate  of  atmospheric  access  to  the 
stores  of  available  energy  and  water  is  also  important  Commonly,  the  mechanisms 
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Most  of  the  many  SVATS  now  available  for  use  in  GCMs  include  calculations  which 
involve  these  key  characteristics,  although  some  (Shuttleworth,  1992)  still  follow  the 
approach  of  many  practicing  hydrologists  in  not  including  explicit  description  of  the 
stomatal  control  and  rainfall  interception  mechanisms  when  calculating  evaporation  rate. 


currently  available  have 
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using  a set  of  simple  aggregation  rules  for  each  parameter  (Shuttleworth,  1991,  1992). 
The  philosophy  is  that,  although  such  simple  aggregation  rules  may  not  stand  the  test  of 
time,  once  in  place  any  critical  evaluation  of  them  based  on  observation  or  modeling  can 
only  result  in  their  replacement  by  superior  alternatives,  and  progress  will  automatically 
follow. 


SVATS  are  required  to ‘describe  area-average  surface  fluxes,  primarily  of  energy  and 
water,  although  some  attempt  is  now  being  made  to  include  description  of  the  exchange 


Let  us  take  the  example  of  momentum  flux,  r , and  the  controlling  parameter,  z#  tne 
aerodynamic  roughness.  In  this  case,  the  relevant  relationship  between  momentum  flux 
and  aerodynamic  resistance  (neglecting  other  surface  parameters  ) is: 


where  zj,  is  a hypothetical  'blending  height'  at  some  level  in  the  air  above  the  ground 
where  turbulent  mining  is  sufficient  for  it  to  be  assumed  that  the  atmosphere  has,  in 
effect,  lost  its  knowledge  of  the  differing  types  of  land  cover  on  the  ground  below. 
Mason  (1988)  estimated  z*  to  be  of  the  order  L/100,  where  L is  the  characteristic 


change  in  temperature;  bc  is  tne  amount  oi  precxpuauon  wmcn  cau  uc  suncu  «u  me 
surface  of  vegetation  (and  soil)  during  rain;  K*  is  the  thermal  conductivity  of  the  soil;  C, 


Table  1.  ’Strawman’  Aggregation  Rules  for  Key  Parameters  in  Soil-Vegetation- 
Atmosphere  Transfer  Schemes  (taken  from  Shuttleworth,  1991;  1993) 


is  the  specific  heat  of  the  soil;  and  Ss  is  the  amount  of  water  which  can  be  stored  per 
unit  area  within  the  vegetation’s  rooting  zone. 
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dependence  of  this  parameter.  The  dependency  of  g,  on  other  variables  was  assumed 

« 

to  be  the  same  or,  in  the  case  of  soil  moisture,  measured.  This  study  gives  similar  good 
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where  f(  ) indicates  a ’stress  function’  for  vapor  pressure  deficit,  de,  temperature,  T,  and 
soil  moisture  deficit,  j,  respectively;  F0  is  the  incident  flux  of  Photosynthetically  Active 


Radiation  (PAR),  while  wF  is  a factor  relating  g,  to  F0.  Analysis  of  the  FIFE  data,  which 
involves  first  removing  the  effect  of  the  stress  functions  to  determine  vF,  reveals  that  the 
value  of  this  last  factor  bears  an  approximately  linear  relationship  to  the  Simple  Ratio,  , 
SR,  with  the  general  form: 
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lative  importance  they  a priori  must  assign  to  the  continental  precipitation 
aporation  fields  and  the  inherently  different  natural  geographical  elements  on  * 


For  small  watersheds  with  dimensions  on  the  order  of  10  km,  within-catchment  recycling 
of  water  is  a small  component  in  the  water  budget  and,  at  this  scale,  the  input  required 
for  hydrological  modeling  is  the  net  input  of  water  across  the  small  catchment.  Given 
such  an  input  for  each  individual  small  basin  in  a larger  watershed,  conventional  and 
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represent  as  a vector  dA.  which  is  defined  to  be  normal  to  the  area  always  pointing  in 
the  outward  direction  from  the  watershed,  as  in  Figure  4c  and  Figure  4d.  The 


* 


storage  components  in  the  watershed. 


The  above  logic,  the  Gauss  theorem,  can  be  used  to  relate  flux  divergence  within  a 
control  volume  to  the  integral  of  flux  over  the  control  surface.  However,  if  we  wish  to 
examine  the  time  rate  of  change  of  a fluid  passing  through  a control  volume,  it  is 


where  the  term  dB/dt  is  the  material  derivative  of  the  fluid  property  (that  is,  its  total 
. rate  of  change  with  respect  to  time);  while  the  first  term  on  the  righthand  side  is  the 


time  rate  of  change  of  the  amount  of  property  stored  in  the  volume. 


eorem. 


Equation  (9)  can  be  accounted  for  by  interpolating  the  vertically  integrated  values  of 
specific  humidity  into  a subgrid  across  the  watershed,  areally  integrating  these  values  at 
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The  time  sequence  of  evolving  global  meteorological  fields  generated  by  the  means  just 
outlined  is  sometimes  referred  to  as  four-dunensional  (4~D)  model  assimilation  data  and 
is  routinely  stored  by  the  major  numerical  weather  prediction  centers  several  times  each 
day  as  a byproduct  of  their  operation.  As  stated  above,  although  substantially 
determined  by  the  assimilated  input  data,  they  are  also  significantly  influenced  by  the 
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data  a«d  the  catchment  scale  to  which  it  is  applied;  or  . 
shortcoming  in  the  GCMs  physics  associated  with  poor  recognition 
of  topography-generated  vapor  extraction  arising  from  processes 
which  occur  at  subgrid  scales. 
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context  of  GCIP,  an  investigation  could  be  made  to  establish  the  scaling  characteristics  ' 
of  the  RTT-based  interface.  This  would  involve  applying  it  to  progressively  larger 
catchments  in  the  Mississippi  River  basin,  each  being  comprised  of  smaller  watersheds 
to  which  the  method  can  also  be  applied  and  for  which  runoff  data  are  available.  It  is 
to  be  anticipated  that  there  will  be  substantial^’ noise’  in  the  relationship  between  RTT 


calculations  and  obser^Qd  runoff  at  the  small  basin  scale  even  in  monthly  average  values, 
but  that  this  will  progressively  reduce  and  then  stabilize  as  the  dimensions  of  the 
sampled  watershed  approach  and  then  exceed  those  of  the  GCM  grid. 
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most  commonly  used  SVATS,  the  Biosphere-Atmosphere  Transfer  Scheme  ( BATS ; Dickinson  at  aL 
1993)  is  evaluated  using  the  surface  energy  and  water  fluxes  collected  in  this  study. 
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After  several  weeks  without  predpitatioa,  on  September 23, 1993,  a bulk  soil  density  of  1700 


\V\ 


to  the  latent  heat  flux,  is  assumed  to  be  proportional  to  these  differences,  thus: 


Ly  exacting  sec  of  criteria  was  applied  to  select  among  these  Bowen  ratio  data  to 


The  diurnal  pattern  of  eddy  covariance  energy  fluxes  are  shown  in  Figure  6(a)  for  a cypic 
y day  (July  22, 1993)  with  clear  sides,  and  in  Figure  6(b)  for  a day  with  substantially  more  cioi 
>ver  and  a late  afternoon  storm  (August  7, 1993).  Typically,  during  dry  periods,  sensible  beat  ns 


In  this  study,  the  initial  ttkkM  validation  run  parameters  appropriate  for  the 
southwestern  USA  within  the  National  Center  for  Atmospheric  Research’s  ( NCAR ) Community 
Climate  Model  version  2 (CGJtfZ),  specifically  standard  parameters  for  vegetation  type  11,  Le.* 


In.  the  model  calibration,  a stand-alone  version  of  the  BATS  was  ran  with  two  secs  of 
parameters.  The  first  set  consisted  of  the  standard  vegetation  and  soil-related  parameters  assigned 
to  the  BATS  when  it  was  used  to  describe  the  semi-arid  southwestern  USA  in  the  COM2,  hi  the 


'O 


Ia  the  12  BATS  soil-rclaicd.  parameter  sets,  the  values  of  and  3 are  explicitly  tied  to  the 


values  of  ocher  soil-related  parameters  through  the  assumpdoa  chat  there  is  a maximum  soil  water 
tension  (of  15  mbar)  against  which  plants  can  obtain  water.  Plant  wilting  behavior  is  thus  formulated 
as  a soil  property,  with  a fixed  relation  to  other  soil  properties  (such  as  hydraulic  conductivity);  and 
in  this  way,  a sec  of  soil  classes  is  determined  by  the  porosity  of  the  sod.  On  this  basis,  the  soil  at 
the  field  site  was  assigned  to  be  the  BATS  soil  texture  class  9,  as  described  above. 


In  comparison  with  the  Bowen  ratio  observations,  the  latent  heat  flux  is  generally 
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enium,  sensible  heat,  water  vapour  and  carbon  dioxide.  Submitted  to  Journal  of 


Webb,  E.  K_  G.  L.  Pearman,  and  R.  T enmng,  1980:  Correction,  of  flux,  measurements  for  densii 
effects  due  to  heat  and  water  vapor  transfer.  Quart.  J.  R.  Meteorological  Soc w 106, 85-10* 


Table  1.  Meteorological  conditions  at  the  study  site  for  the  typical  dry  and  wet  days  shown  in 
Figures  6(a),  6(b),  7(b)  and  7(c). 


Max.  Wind  Speed  (m  s'*)  4_5  9.9 

Total  Daily- Precipitation  0.0  17 _5 

(mm) 


Table  3.  The  BATS  soil-related  parameters  as  used  in  the  National  Center  for  Atmospheric 
Research  Community  Climate  Model  version  2 (COM2)  for  the  “semi-dizerf1  areas  of 
southwest  USA  ( BATS  soil  texture  class  3)  and  the  modified  site-specific  parameters 
which  produced  an  improved  model  description  of  the  observed  data. 
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Research  Note 


Abstract 


used  as  forcing  variables  for  an  equivalent  array  ot  stand-alone  mospnere-Aunospnere  iransier 
Scheme  (BATS)  models  to  describe  the  evolution  of  soil  moisture  and  surface  energy  fluxes  in 
response  to  the  prevalent,  heterogeneous  pattern  of  convective  precipitation.  The  calculated  area- 
average  behavior  was  compared  with  that  given  by  a single  aggregate  BATS  simulation  forced  with 


rea-average  meteorological  data.  Heterogeneous  rainfall  gives  rise  to  significant  but  partly 
ompensating  differences  in  the  transpiration  and  the  intercepted  rainfall  components  of  total 
vaporation  during  rain  storms.  However,  the  calculated  area-average  surface  energy  fluxes  given 


This  research  note  addresses  the  question,  “In  the  extreme  case  of  the  semi-arid  environment 


landscape,  the  performance  of  the  hypothetical  aggregation  rules  was  compromised.  This  result  is 
explained  by  the  non-linear  response  of  vegetation  to  changing  soil  moisture  in  the  BATS  and,  it  is 
assumed,  in  nature  as  well  — see  Arain  et  al.  (1996)  for  greater  detail. 

An  analysis  of  satellite-derived  land  cover  classes  shows  that  the  proportion  of  irrigated  land 


Draft  (1.6) 

falling  within  an  individual  grid  square  used  in  meteorological  models  is  rarely  greater  than  10%  in 


aggregation  rule  performance  exists  in  natural  situations.  This  note  explores  this  specific  issue  for 
the  deliberately  extreme  case  of  soil  moisture  heterogeneity  resulting  from  intermittent  localized 
convective  rain  falling  in  a semi-arid  region  of  the  southwestern  USA. 


2.  Field  Data 

The  data  used  in  this  study  is  from  the  Walnut  Gulch  Experimental  Watershed  (31°43’N, 
1 10°00W),  which  is  located  120  km  southeast  of  Tucson,  Arizona,  and  which  comprises  the  upper 


convective  precipitation  during  a summer  monsoon  season  encompassing  the  months  of  July  through 


Draft  (1.6) 


with  1 1 runoff  measuring  flumes,  and  two  sets  of  weather  and  energy  flux  measuring  stations 


2.1.  Synthesis  of  Distributed  Data  Arrays 


Draft  (1.6) 

maximum.  Some  larger,  short-term  differences  were  observed  (particularly  in  temperature  and  vapor 


In  consequence  of  the  just-descnbed  observations  ana  moaenng  resuus,  u is  appropriate  to 
assume  that  the  values  of  incoming  solar  radiation,  near-surface  air  temperature  and  vapor  pressure 
deficit  required  by  the  stand-alone  version  of  the  BATS  are,  in  effect,  constant  .across  the  26.9  by 
12.5  km  study  area.  For  reasons  of  data  consistency,  the  preferred  values  of  meteorological  forcing 


rain  gauge  network  (Haider,  1994)  have  demonstrated  that  the  multiquadric-biharmonic  technique 


Draft  (1.6) 

is  the  preferred  method  of  interpolating  precipitation  because  it  combined  practicality  with 


study  site. 

Performance  of  these  two  parameter  sets  was  evaluated  by  comparing  the  resulting  modeled 
surface  fluxes  with  the  measurements  at  the  two  observation  sites.  However,  the  value  of  this 


comparison  was  restricted  by  known  weaknesses  in  the  flux  measuring  instruments  used  at  these 
sites.  The  observation  of  sensible  heat  flux  was  made  at  a height  of  9 m using  the  eddy  correlation 
method  with  a single  “Gill”  propeller  to  measure  the  vertical  wind  velocity.  A previous  comparison 


2 (from  the  8 available  color  classes  ranging  from 
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for  definition  of  these  classes.  Simulations  using  these  parameters  provided  poor  agreement  with 


minimum  leaf  area  indices),  and  many  soil  parameters  (including  porosity,  saturated  hydraulic 
conductivity,  and  wet  and  dry  soil  albedos).  In  addition,  a prior  modeling  study  earned  out  at  a 
Sonoran  Desert  site  near  Tucson,  Arizona  (Unland  et  al.,  1996)  suggested  that  two  other  parameters 


for  latent,  sensible,  and  ground  heat  flux  respectively. 


. Modeling  Studies 


vegetation  temperatures  for  the  two  most  extreme  elements  in  the  model  array  were  selected  and 


forcing  variables  that  had  been  modeled  to  occur  at  this  height  in  the  BATS-ABL  run).  In  this  way 


an  adequate  area-average  description.  In  this  way  the  area-average  interaction  with  the  realistic 
representation  of  precipitation-induced  heterogeneity  in  soil  moisture  can  be  compared  with  that 
given  when  heterogeneity  is  neglected  in  the  modeling  process.  If  the  two  outputs  disagree  with  one 
another,  then  failure  of  aggregation  rules  (analogous  to  that  which  Arain  et  al.  (1996)  reported  as 
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equivalent  runs  using  the  site-specific  parameters,  none  of  the  conclusions  given  below  are  affected 
because  these  conclusions  are  based  on  the  difference  between  aggregate  and  distributed  array  runs 


Figure  2 shows  a comparison  between  the  calculated  latent  heat,  sensible  heat  and  ground 
heat  fluxes  given  by  the  aggregate  BATS  model  and  those  given  by  the  distributed  array  of  BATS 


vaporation  flux  for  the  study  area  that  originate  (a)  as  transpiration,  and  (b)  as  the  evaporation  of 


On  the  basis  of  the  above  comments,  it  is  clear  that  heterogeneous  rainfall  can  result  in  a poor 
description  of  short-term  evaporation  rates  during  and  immediately  after  rainfall  because  the 
description  of  the  rainfall  interception  process  is  imperfect  in  a one-dimensional,  aggregate 
representation.  This  feature  has  been  recognized  for  some  time,  see,  for  example,  Shuttleworth 
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(1988).  However,  in  the  case  of  semi-arid  regions,  where  there  is  little  runoff  and  where  water 


The  basic  reason  for  this  is  that  the  amount  of  rain  in  an  individual  precipitation  event  is 


at  any  time,  only  about  10%  of  the  non-linear  relationship  between  soil  moisture  am 
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evaporation  is  being  sampled  within  a given  soil  moisture  pattern,  and  the  non-linear  relationship 


Arain,  A.M.,  1994.  Spatial  aggregation  of  vegetation  parameters  in  a coupled  land  surface- 
atmosphere  model.  M.S.  Thesis,  The  University  of  Arizona,  Tucson,  USA.  150  pp. 

Arain,  A.M.,  Michaud  J„  Shuttleworth,  W.  J„  and  Dolman,  A.  J.,  1996.  Testing  of  vegetation 
parameter  aggregation  rales  applicable  to  the  Biosphere- Atmosphere  Transfer  Scheme  (BATS)  and 
the  FIFE  site.  J.  Hydrol.,  177:1-22. 
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Figure  Captions 


area-average  meteorological  data  (full  line)  and  the  distributed  array  of  BATS 
models  forced  with  distributed  data  (broken  line)  and  (c)  the  interpolated 
precipitation  from  day  of  year  213  to  215  of  1990  (August  1-  August  3, 
1990).  There  were  7 localized  rain  storms  during  this  period. 
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Table  1 Modeled  Water  Balance  Components  for  the  Aggregate  and 

Distributed  BATS  models  using  site-specific  parameters. 
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Abstract 


heat  and  water  vapour  in  heterogeneous  terrain.  It  is  suggested  that  this 
simple  parameterization  can  adequately  take  into  account  the  effect  of 
variation  in  surface  cover  on  the  fluxes  of  heat  and  water  vapour. 
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surface  cover  on  the  fluxes  of  heat,  momentum  and  water  vapour,  but  would 
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neglect  advective  processes.  This  paper  is  primarily  concerned  with  the 
theoretical  approach,  which  tries  to  improve  on  the  previously  mentioned 
studies  by  taking  advective  processes  into  account. 
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parameterizations  based  on  the  concept  of  sparse  crop  energy  balance 
models  (Dolman,  1993a,  Blyth  and  Dolman,  1994)  are  then  used  to  show  the 


the  frictional  velocity,  u.,  the  appropriate  velocity  scale  at  canopy  or  surface 
layer  level.  Raupach  (1992)  presents  a framework  to  estimate  the  relevant 
horizontal  length  scales  involved  in  this: 
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It  is  possible  to  extent  this  sparse  crop  framework  to  larger  areas,  by  following 
the  same  strategy  of  applying  a one-dimensional  framework  to  solve  for  a 


and  heat  are  still  approximately  in  local  equilibrium,  whereas  the  mean  proper- 
ties such  as  windspeed  and  concentration  are  sufficiently  well  mixed  to 
assume  that  they  are  (horizontally)  homogeneous  (Wieringa,  1986,  Mason, 
1988).  Thus,  for  small  scale  land  cover  variability  the  relevant  scales  are  the 
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the  fluxes  are  still  approximately  in  local  equilibrium  whereas  the  windspeed 
represents  almost  an  area  average  and  has  "blended  out . Several 


where  L*  is  a characteristic  length  scale  for  the  horizontal  variations.  Several 
other  estimates  of  the  blending  height  may  be  used  for  practical  purposes, 
such  as  the  diffusion  height  scale  or  1^200  (Mason,  1988)  or  the  height  of  the 
lowest  model  layer,  but  equation  4 provides  arguably  the  best  physical 
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<T>  = 


heat  flux,  XE  the  latent  heat  flux,  r a resistance  with  subscript  ah  and  ae  for 
heat  and  water  vapour  respectively,  p the  density  of  air,  cp  the  specific  heat 
of  air,  X the  latent  heat  of  vaporization,  y the  psychrometer  constant  and  T 
temperature,  e water  vapour  pressure  where  the  subscript  s denotes  a surface 


value  and  a a value  at  a reference  level  (blending  height)  in  the  air. 


QO 
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which  show  that  the  effective  resistance  in  heterogeneous  terrain  is  an  area- 
average  weighted  by  flux  (e.g.  Blyth  et  al.,  1993).  This  makes  averaging 
resistances  a non  trivial  task  as  a priori  knowledge  of  the  fluxes  is  required.  It 


is  also  of  interest  to  note  that  the  aerodynamic  and  surface  resistance  for 
latent  heat  appear  as  a single  term  in  the  equations.  It  is  possible  to  extend 
this  analysis  to  the  case  of  flux  equations  with  roughness  lengths  (Wood  and 
Mason,  1991),  but  the  conclusions  remain  essentially  the  same. 
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conclusions  were  obtained  by  Wood  and  Mason  (1991)  and  Blyth  et  at.  (1993) 
and  Dolman  (1993b). 


These  numerical  studies  have  been  useful  to  look  in  detail  at  some  of  the 
processes  relevant  for  aggregation.  By  dividing  up  the  flow  in  a model 


(0 


The  main  conclusions  from  the  numerical  models  studies  are  that  the  velocity 
perturbations  are  generally  of  greater  magnitude  than  the  scalar  perturbations 
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certain  height,  the  blending  height.  They  then  used  various  resolutions  of  the 
model  to  test  whether  simple  averaging  rules  (Dolman,  1993a)  could  provide 


use  simplified  representations,  effective  parameters,  in  mesoscaie  moaeis  as 
long  as  the  original  surface  information  is  sufficient.  This  conclusion  is  further 
based  on  the  fact  that  the  advective  processes  operating  in  such  a 
heterogeneous  area  cancel  each  other  to  a sufficient  extent  to  have  a 


negligible  effect  on  the  area  average  fluxes. 
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PARAMETERIZATIONS 
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^ V 


where  k is  van  Karman's  constant  and  the  integral  stability  correction  for 
momentum.  At  the  blending  height,  the  horizontal  velocity  is,  by  definition, 
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In  contrast  to  momentum  transport,  transport  of  heat  lacks  a no  slip  condition 


J1 

H-8 


Averaging  procedures  for  heat  and  momentum  are  therefore  not  necessarily 
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(Monteith,  1981)  which  consist  of  a combination  of  the  surface  energy  balance 
and  the  transfer  equations.  For  latent  heat: 
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Table  1 shows  the  results  of  this  comparison  for  surface  resistance  whilst 
Table  2 shows  results  for  aerodynamic  resistances  to  sensible  and  latent  heat 
respectively  for  the  blending  height  averaging  only.  If  the  range  in  surface 


CD 


Han  Dolman,  Version  02,  5 July  1995 


THE  ROUGHNESS  LENGTHS  FOR  HEAT  AND  MOMENTUM 
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elegantly  by  parameterizing  the  laminar  boundary  layer  resistances,  which  are 
felt  only  by  scalar  transport,  explicitly  (Owen  and  Thompson,  1963). 
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better  to  use  a surface  energy  balance  model  which  explicitly  models  the 
different  sources  of  heat  momentum  and  water  vapour.  A similar  conclusion 
may  be  drawn  for  areas  with  large  heterogeneity. 


DISCUSSION 


2-D  numerical  boundary  layer  model.  The  agreement  is  surprisingly  good, 
suggesting  that  the  concept  of  blending  height,  coupled  with  adequate  descrip- 
tions of  the  land  surface,  can  reliably  estimate  surface  fluxes  from  heterogen- 
eous surface  cover  at  a range  of  length  scales. 


In  principle  by  using  the  (sparse  crop)  surface  energy  balance  models  at 
larger  scales  than  for  which  they  were  originally  developed-as  shown  in  the 
previous  section  of  this  paper-  the  need  for  deriving  effective  values  vanishes 
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Simple  averaging  of  the  series  and  parallel  resistances  is  a good  option  to 
estimate  effective  parameters  for  resistances.  When  more  detail  is  required  the 
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Question  1 , dr.  A.  Jochum,  DLR,  Germany. 
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Question  2 Dr.  M.R.  Raupach,  CSIRO,  Australia. 


remote  sensing  images  may  be  the  way  forward  to  determine  the 
shape  of  these  spectra.  At  this  stage  the  most  probable  solution  to  the 
problem  is  to  classify  the  length  scales  along  the  lines  suggested  in  this 
paper.  A natural  upper  limit  exists  at,  say,  10  km,  as  from  there  on  the 
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where  VE  is  evaporation,  G soil  heat  flux,  H sensible  heat  flux 
and  the  subscripts  refer  to  the  resolution  of  the  land  surface 
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Figure  8.  Ratio  of  resistance  from  blending  height  to  model  layer  to  that 
of  canopy  level  to  blending  height  (see  Figure  5 and  text  for 
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The  scaling  characteristics  of  soil  parameters: 

From  plot  scale  heterogeneity  to  subgrid  parameterization 
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between  soil  textural  composition  and  hydraulic  characteristics,  these  aggregatod 
parameters  predict  evaporation  fluxes  well,  but  fail  to  predict  water  balance  terms 
such  as  percolation  and  runoff.  This  is  a serious  drawback  which  could  eventually 
hamper  the  improvement  of  the  representation  of  the  hydrological  cycle  in 
mesoscale  atmospheric  models  and  in  GCMs. 
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with  Z the  centre  of  the  soil  profile  depth.  The  area  average  values  of  the  soil 
moisture  content  8 and  of  the  hydraulic  conductivity  K are  dependent  on  the  local- 
scaie  solution  of  the  soil  suction  h (x,t) : 

(e)(X,0  = (0{/j(x,f))>  (7 

(K){X,f)  = (K{h{x,1)))  (8 
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3.  Averaging  of  saturated  hydraulic  conductivity 

The  saturated  vertical  hydraulic  conductivity  K A varies  strongly  between  soil  types 
(Clapp  and  Homberger,  1978;  Cosby  et  al,  1984,  Wosten  et  al,  1995).  In  predictive 
models  of  soil  hydraulic  conductivity  (Van  Genuchten  et  al,  1 992;  Clapp  and 
Hornberger,  1978  and  Mualem.1976,  among  others)  Ka  is  generally  treated  as  a 


measurable  parameter.  Ks  determines  to  a large  extend  the  partitioning  of  rainfall 
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1984).  In  physical  systems,  three  types  of  similarity  can  be  defined:  geometric. 


kinematic  and  dynamic.  In  this  paper,  geometric  and  dynamic  similarities  will  be 
discussed. 


x and  of  a reference  soil.  From  viscous  flow  theory  (Raats,  1990)  the  following 
scaling  rules  can  be  derived  for  the  soil  moisture  content  0,  pressure  head  h and 
hydraulic  conductivity  K : 


These  equations  imply  that  the  soil  water  retention  and  hydraulic  conductivity 
characteristics  at  a given  water  content  9X  at  any  location  can  be  related  to  mean  h 
- 0 1 - K functions,  and  are  of  a general  form  : 


Values  of  exponent  n in  scaling  equation  (18)  for  several  variables  and 
parameters  used  to  model  soil  water  flow.  D is  the  water  diffusivity  ( D = k dh/db), 
MFP  is  a matrix  flux  potential  (Shaykewich  and  Stroosnijder,  1977)  and  vis  the 
velocity;  for  other  variables  see  text. 
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variables  can  be  converted  into  dimensionless  form  (Kuhnel,  1989): 


two  universal  scale  factors  can  be  derived  for  both  time  and  depth  dependent  soil 
water  flow  : 
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There  are  several  closed  analytical  solution  for  water  diffusivity,  such  as  (Fujita 
1952): 


5.  Effective  soil  parameters  from  inverse  modelling 
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Volumetric  water  content  of  the  profiles  on  selected  days; 
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6.  Aggregated  soils 

Techniques  to  estimate  representative,  aggregated  soil  hydraulic  characteristics  of 
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area  equal  to  4/9,  sand  having  the  fractional  area  of  3/9  and  the  loam  with 
fractional  area  equal  to  2/9.  Soil  hydraulic  characteristics  are  described  by  the 
Clapp  and  Homberger  (CH)  (1978)  model: 
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To  each  of  the  soil  types,  mean  clay  and  sand  percent  compositions  ( C,  and  SJ 
have  been  assigned  according  to  Clapp  and  Homberger  (1978).  For  example,  the 
textural  soil  type  "day”  has  a mean  clay  fraction  equal  to  0.63,  and  a sand  fraction 
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The  techniques  presented  in  this  paper  to  scale  and  aggregate  the  soil 
characteristics  all  potentially  qualify  for  use  in  large  scale  meteorological  models. 
One  of  the  most  attractive  and  perhaps  surprising  results  is  the  effective  behaviour 
of  the  reference  curve,  obtained  through  similar  media  scaling  (section  4).  Further 
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Figure  10.  Results  of  the  SWAP/MITRE  simulations  with  aggregated  and 
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set  of  functions  which  depend  mainly  on 


depend  on  the  flow,  or  whether  one-dimensional  or  three-dimensional 
representations  are  used. 


Q4  (dr.  Detlef-Schultze): 

In  averaging  soil  water  you  cannot  neglect  plant  cover,  because  it  will  determine 
porosity  and  the  amount  of  deep  drainage.  Grasses  and  shrubs  have  different  root 
structures  and  determine  soil  properties  independent  of  soil  types, .like  clay,  sand 
and  loam. 


Distance  (m) 


log  (-h)  [cm] 


log  (-ft)  (cm) 


% Sand 


Reference 
US  Profile 
■I  Evaporation 


240 


Mean 

Reference 


Time  (days)  1982 


■ Sand 


log  K (cm/day) 


With  increasing  interest  in  airborne  and  satellite-based  sensors  for  mapping  regional 


and  global  energy  balance,  there  is  a need  to  determine  the  uncertainty  involved  in 
aggregating  remotely-sensed  variables  [surface  temperature  (T^)  and  reflectance  (p)]  and 


pixel  resolution  and  averaged  to  the  coarser  resolution;  and  second,  calculated  directly  at  the 


coarse  resolution  by  aggregating  the  fine-resolution  data  to  the  coarse  scale.  Results  showed 
that  the  error  in  the  aggregation  of  Tk  and  p was  negligible  for  a wide  range  of  conditions. 

4 

However,  the  error  in  aggregation  of  H and  AE  was  highly  influenced  by  the  heterogeneity 
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INTRODUCTION 


approaches  require  some  supplementary  meteorological  and/or  surface  information. 

The  extent  of  surface  information  required  draws  a distinction  between  those 
approaches  that  can  be  applied  at  a local  scale  and  those  limited  to  regional  application 


(Moran  and  Jackson,  1991).  Most  local  scale  methods  rely  on  site-specific  measurements 
of  aerodynamic  and  atmospheric  conditions  and  apply  only  to  an  area  over  which  the 
ground-based  measurements  can  be  extrapolated.  Regional  scale  methods  generally  rely  on 


o 


(particularly,  radiometric  temperature,  Tr)  and  aggregation  of  energy  balance  components 


pixel-resolution  values  of  radiance  were  averaged  to  the  coarser  scale  and  then  a value  of  H 
was  computed.  The  difference  between  these  two  computations  of  H revealed  the  error  due 
to  non-linearities  in  the  relation  between  surface  radiance  and  H. 


The  presentation  is  organized  into  sections,  with  the  first  section  presenting  the  basic 
energy  balance  theory  that  is  the  foundation  for  most  local-  and  regional-scale  models  and 
defining  the  relation  between  kinetic,  radiometric  and  aerodynamic  temperatures.  The 


where  R^,  G,  H and  XE  are  in  units  of  W m'2,  and  G,  H and  XE  are  positive  when  directed 
away  from  the  surface.  A.E  is  a function  of  evaporation  rate,  E (kg  s m ),  and  heat  of 
vaporization,  L (J  kg'1)),  but  is  typically  found  as  a residual  in  Eq.  [1]. 


computed  using  remotely-sensed  measurements  to  estimate  a and  RLI,  and  ground-based 
instrumentation  to  measure  the  incoming  terms  (Jackson,  1984;  Brest  and  Goward,  1987). 
The  G term  is  dependent  upon  the  gradient  of  temperature,  T (°K),  with  soil  depth. 
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where  Cv  the  volumetric  heat  capacity  of  air  (J  °C' m*3),  and  ra  is  the  aerodynamic  resistance 
to  heat  transport  (s  m'1).  This  resistance  for  neutral  conditions  (where  T0  = TJ  is 


where  z0  and  d are  roughness  length  and  zero-plane  displacement  (m),  respectively,  U is 
wind  speed  (m  s'1),  z is  the  height  (m)  above  the  surface  where  u is  measured,  and  k (=0.41) 


neglecting  atmospheric  effects  in  measurements  acquired  with  an  aircraft  or  satellite-based 
sensor)  in  a finite  spectral  band  is  related  to  the  emittance  of  the  surface  (Mk,  W m ) with 
the  expression 


= eMk  + (l-e)Ms 


here,  the  Stefan-Boltzman  expression  is  substituted  for  each  term  in  Equation  6 to  obtain  the 


In  the  above  expression,  the  term  B includes  the  incoming  sky  radiation  in  the  spectral  band 
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- a:  There  is  a great  deal  of  published  evidence  that  the  relation  between  radiance  and 
reflectance  (and  a)  is  linear  for  most  sensors  and  most  atmospheric  conditions  (Hall 


Again,  since  previously-published  studies  reported  a linear  relation  between 
reflectance  and  radiance,  there  is  little  chance  for  error  in  the  aggregation  of  values 
of  G from  local  to  regional  scales. 


- XE:  Since  XE  is  computed  as  a residual  in  the  energy  balance  equation  (Eq.  [1]),  the 


The  relation  between  radiance  measured  by  a satellite  sensor  and  values  of  surface 


temperature  is  not  necessarily  linear.  For  example,  the  relation  between  Tr  and  spectral 
radiance  (L)  is  given  for  Landsat  Thematic  Mapper  (TM)  as 


derived  from  radiance  observed  with  sensors  of  low-  and  high-resolution  over  the  same  area 


1 .0094+0 .042(ln(ND VI)) , 


where  NDVI  is  the  Normalized  Difference  Vegetation  Index  [(PNiR-Pred)/ (PNiR+Pred)]  and  Pnir 
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Third,  there  is  a non-linear  relation  between  Tk  and  T0  over  a range  of  vegetation 
densities.  That  is,  the  difference  between  T0  and  Tk  is  small  for  fully-vegetated  surfaces  (or 
sites  with  little  difference  between  soil  and  vegetation  temperatures)  and  increases  rapidly 


with  increasing  heterogeneity  of  the  surface  (when  the  difference  between  soil  and 


(Brutsaert,  1975b).  This  results  in  additional  resistance  to  heat  transfer  that  can  be  expressed 
in  the  resistance  equation  by  including  different  roughness  lengths  for  momentum  (subscript 


included  in  the  resistance  computation,  resulting  in  a non-linear  relation  between  (T^TJ  and 


where  s^  is  site-specific  and  was  computed  as  0.17  for  Owens  Valley  rangeland  and  0.13  for 


Walnut  Gulch  rangeland  (Kustas  et  al.,  1994b). 


during  the  "dry"  season  than  during  the  "wet"  season  when  vegetation  and  soil  moisture  are 
highly  variable.  However,  an  early-morning  measurement  during  the  wet  season  would  have 
less  variability  than  a mid-day  measurement  during  the  dry  season  (due  to  differential 


heating  of  vegetated  and  bare  surfaces).  Furthermore,  a fine-resolution  image  during  the  dry 


theoretical  background  which  will  be  provided  here,  though  readers  are  encouraged  to  refer 
to  previously-published  works  (Price,  1990;  Gillies  and  Carlson,  1994;  Nemani  and  Running, 


A scattergram  of  the  spectral  vegetation  index  (e.g.  NDVI  or  soil-adjusted  vegetation 
index  (SAVI),  Huete,  1988)  vs.  Tk-Ta  generally  results  in  a trapezoidal  shape  that  is 
indicative  of  the  variation  in  vegetation  cover  and  evapotranspiration  (e.g.,  Figure  1). 


possible  surface  temperatures  for  a range  of  vegetation  cover  from  0 to  100%  using  a 


Vegetation-Index/Temperature  (VIT)  Trapezoid  that  encompassed  the  maximum  possible 
variability  in  the  vegetation  index  and  values  of  (T^-TJ  for  one  site  on  one  date  (Figure  1). 
The  VIT  Trapezoid  theory  defines  Tk-Ta  for  four  extreme  situations  by  changing  the 


1)  lull-cover,  well-watered  vegetation,  where  rc  is  the  canopy  resistance  at  potential 


o 


The  left  edge  of  the  trapezoid  was  related  to  conditions  of  potential  evapotranspiration 


In  order  to  test  any  schemes  for  aggregation  of  Tr  and  H,  variability  in  the  following 


Most  of  the  analysis  was  limited  to  eight  sites  within  the  watershed,  termed 
METFLUX  stations  (herein  referred  to  as  MF1  to  MF8),  containing  instrumentation  for 
measuring  both  general  meteorological  conditions  and  estimating  the  surface  energy  balance 


(Kustas  et  al.,  1994a).  The  METFLUX  sites  were  located  along  two  parallel  transects 


resulting  in  very  heterogeneous  conditions. 


delimited  over  which  surface  temperature  and  reflectance  were  measured  from  a 


(5 


conditions.  Thermal  data  were  selected  for  Kendall  and  Lucky  Hills  on  three  dates 
at  0.09  km  AGL:  DOY  209  @ 10:12,  DOY  209  @ 14:40  and  DOY  221  @ 7:30. 
Thermal  data  acquired  at  higher  altitudes,  0.92  km  and  3.35  km  AGL,  were  selected 


for  Lucky  Hills  on  DOY  222  @ 14:11  and  Kendall  on  DOY  209  @ 14:51, 


approximately  2.3  by  2.9  km  were  extracted  for  areas  surrounding  each  of  the  8 
METFLUX  sites.  Though  the  NS001  sensor  provides  data  in  8 spectral  bands,  at  this 
time  only  the  thermal  data  were  available.  We  limited  our  analysis  to  data  collected 


m 

(N 


We  emphasize  that  these  data  sets  were  selected  for  analysis  because  they  covered  a 


is  at  the  regional  and  global  scales,  the  aggregation  of  our  finest  resolution  data  (less  than 
plant  diameter)  provided  the  opportunity  for  insight  into  the  uncertainty  of  aggregation  at 
coarser  scales  (e.g.,  sites  characterized  by  distinctive  patch  vegetation). 


f T and  H were  each  computed  in 


pixel-resolution  values  of  radiance  were  averaged  to  the  coarser  scale  (e.g.  1 km  for 


pixel-resolution  values  of  Tk  were  averaged  to  the  coarser  scale,  and  then  a value  of 
H was  computed  using  Eqs.  [10]  - [12]  based  on  average  values  of  Tk.  All  tests  were 
conducted  using  resistance  equations  that  included  the  kB-1  factor,  using  Sj^O.17. 
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trapezoid,  indicating  higher  XE  values  associated  with  the  wet  season  (Figure  3b).  The 
general  variability  of  NDVI  and  Tk-Ta  was  much  greater  than  for  DOY  156  at  both 
resolutions.  The  variability  of  the  yoke  data  at  Lucky  Hills  (MF1)  was  greater  than  the  same 


data  for  Kendall  (MF5).  Similar  to  results  for  DOY  156,  the  variability  of  the  yoke  data  was 


Hills  and  Kendall  (Figure  4a  and  Table  2).  This  was  due  to  the  differential  thermal 
inertia  and  evapotranspiration  of  vegetation  and  soil,  and  it  generally  results  in  greater 
absolute  variability  near  noon.  However,  since  one  of  the  issues  to  be  addressed  in 


variability  than  the  finer  resolution  data;  the  range  of  Tr-Ta  values  decreased  from  27 
°C  (at  0.2  m resolution)  to  13  °C  (at  6.0  m resolution).  This  was  the  case  even 
though  the  coarser-resolution  image  covered  a ground  surface  area  that  was  larger  and 


more  variable  than  the  smaller  area  covered  by  the  finer-resolution  image. 
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In  previous  sections,  we  identified  several  sources  of  non-linearity  in  the  relation 
between  radiance  and  Tr  (due  to  sensor  calibration),  between  Tr  and  Tk  (due  to  emissivity), 
and  between  Tk-Ta  and  H (due  to  stability  corrections  qf  r and  variability  in  surface 


roughness,  zj.  Here,  we  used  data  from  the  yoke-based  radiometers,  aircraft-based  thermal 


a large  range  of  Tr  values  and  resulted  in  aggregation  errors  less  than  0.05  °C  for  temperature 
ranges  of  30°C  for  prairie  vegetation  in  Kansas  (Humes  and  Sorooshian,  1994)  during  FIFE, 
the  First  ISLSCP  Field  Experiment  (Sellers  et  al.,  1988).  Similar  results  were  found  here  for 


O 


values  using  the  two  aggregation  schemes  was  less  than  0. 1 °C.  The  error  in  aggregation 


eterogeneous  sites,  such  as 


In  the  following  tests,  we  assumed  that  zQ  was  constant  at  both  the  pre-  and  post- 
aggregation scale.  In  this  way,  we  could  investigate  the  variability  in  H associated  with 
variability  in  the  stability  correction  of  vj  and  variability  in  kB  . By  comparing  results  at 


It  is  notable  that  the  site  conditions  during  data  acquisition  of  yoke-based  data  (@ 


For  moderate  resolution  (6.0  m)  NS001  data,  error  in  aggregation  associated  with  the 
stability  correction  was  small  (<  6 W m’2)  for  all  METFLUX  sites.  As  with  the  yoke-based 
data,  the  data  were  acquired  during  mid-day  when  unstable  conditions  dominated  the  site 


(Tk-Ta>0)  for  most  pixels.  Results  were  similar  for  the  coarse  resolution  (120  m)  TM  data. 
The  error  associated  with  aggregation  of  data  acquired  during  the  dry  season  was  negligible 
for  both  the  WGEW  and  RIP  sites.  For  data  acquired  during  the  wet  season  (DOY  252)  at 


fsl 

fO 
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(nearly  50  W m'2)  between  values  of  H aggregated  with  differing  z0  values,  especially  in 
mid-afternoon  (Figure  9).  There  was  negligible  error  (<  5 W m'2)  associated  with 
measurements  in  the  early  morning  when  values  of  H were  small  (-  35  W m 2).  For  NS001 


data  with  moderate  resolution  (6.0  m),  the  error  in  aggregation  of  H was  smaller  (<  25  W m' 


the  use  of  a constant  kB*‘  term  would  minimize  the  errors  in  aggregation  of  H for 


heterogeneous  sites,  it  would  substantially  increase  the  absolute  error  in  estimation  of  H 


rf 


CONCLUSIONS 


This  study  emphasized  the  need  to  account  for  site  heterogeneity  in  selection  of  a 
scheme  for  aggregation  of  surface  energy  balance  components.  This  is  particularly  true  for 
sparsely- vegetated  sites,  such  as  the  Walnut  Gulch  semi-and  rangeland.  However,  the 


quantification  of  site  heterogeneity  is  not  an  easy  task  since  it  can  be  related  not  only  to 
vegetation  cover  but  also  to  time  of  day,  time  of  year  and  sensor  spatial  resolution.  Thus, 
it  is  important  to  understand  the  sources  of  non-linearity  in  the  computation  of  H and  AE  in 


xn 


- Sites  which  are  characterized  by  patch  vegetation,  where  the  pixel  resolution  is  less 


in  surface  roughnesses.  These  results  were  obtained  with  the  aggregation  of  sites  where 
roughness  varied  by  only  0.03  m.  The  aggregation  error  for  sites  with  greater  variation  in 
roughness  lengths,  such  as  mixtures  of  grassland  and  forest,  could  be  many  times  larger  than 


ACKNOWLEDGEMENTS 


<L> 

U 

U 

& 

<u 


Ou 


a 

A 


vO 

ro 


Rob  Perry  who  were  responsible  for  some  aspects  of  data  processing.  We  also  appreciate 


h 

m 


REFERENCES 


England,  p.299. 

Brutsaert,  W.H.  (1986)  Catchment-scale  evaporation  and  the  atmospheric  boundary  layer, 


CO 

i 

GO 

ON 

CO 

(N 


CQ 


oo 


Chehbouni,  A.,  W.D.  Nichols,  J.  Qi,  E.G.  Njoku,  Y.H.  Kerr  and  F.  Cabot  (1994)  On  the  use 


Fuchs,  M.  and  C.B.  Tanner  (1966)  Infrared  thermometry  of  vegetation,  AgronuX  58:597- 


Jackson,  R.D.  (1988)  Surface  temperature  and  the  surface  energy  balance,  IN  Flpw  and 
Transport  in  the  Natural  Environment:  Advances  and  Applications,  W.L.  Steffen  and 
O.T.  Denmead  (Eds.),  Springer-Verlag,  New  York,  pp.  133-153. 


Jackson,  R.D.,  M.S.  Moran,  P.N.  Slater  and  S.F.  Biggar  (1987)  Field  calibration  of  reference 
reflectance  panels,  Rem.  Sens.  Environ.  22:145-158. 

Jacobs,  C.M.J.  and  H.A.R.  de  Bruin  (1992)  The  sensitivity  of  regional  transpiration  to  land- 

4 

surface  characteristics:  significance  of  feedback,  J.  of  Climate  5:683-698. 


Kustas,  W.P.  and  C.S.T.  Daughtry  (1990)  Estimation  of  the  soil  heat  flux/net  radiation  ratio 
from  spectral  data,  Agric.  For.  Meteorol.  49:205-223. 


Kustas,  W.  P.,  D.C.  Goodrich,  M.S.  Moran,  et  al.  (1991)  An  interdisciplinary  field  study  of 


1259. 

Kustas,  W.P,  E.M.  Perry,  P.C.  Doraiswamy  and  M.S.  Moran  (1994c)  Using  satellite  remote 
sensing  to  extrapolate  evapotranspiration  estimates  in  time  and  space  over  a semiand 


rangeland  basin,  Rem.  Sens.  Env.  (in  press). 

Lang,  A.R.G.,  K.G.  McNaughton,  C.  Fazu,  E.F.  Bradley  and  E.  Ohtaki  (1983)  Inequality  of 
eddy  transfer  coefficients  for  vertical  transport  of  sensible  and  latent  heats  during 


Monteith,  J.L.  (1973)  Principles  of  Environmental  Physics,  Elsevier  Press,  New  York,  241 


o 


evapotranspiration  using  remotely  sensed  inputs,  J.  Env.  Oual.  20:725-737. 

Moran,  M.S.,  R.D.  Jackson,  L.H.  Raymond,  L.W.  Gay  and  P.N.  Slater  (1990)  Mapping 
surface  energy  balance  components  by  combining  Landsat  Thematic  Mapper  and 


On 

Tt 

m 


o 

m 

ro 


o 

cn 


<u 

t/5 


0* 


m 

Tf 


Topical  Svmp.  on  Combined  Optical-Microwave  Earth  and  Atmosphere  Sensing.  22- 


Price,  J.C.  (1990)  Using  spatial  context  in  satellite  data  to  infer  regional  scale 
evapotranspiration.  IEEE  Trans.  Geosci.  and  Rem.  Sens.  28:940-948. 

Richardson,  A.J.  and  C.L.  Wiegand  (1977)  Distinguishing  vegetation  from  soil  background 


Sellers,  P.J.,  F.G.  Hall,  G.  Asrar,  D.E.  Strebel  and  R.E.  Murphy  (1988)  The  first  ISLSCP 


arid  areas,  J.  Applied  Meteorol.  (in  press). 

Van  de  Griend,  A. A.  and  M.  Owe  (1993)  On  the  relationship  between  thermal  emissivity 
and  the  normalized  difference  vegetation  index  for  natural  surfaces,  Int.  J.  of  Rem, 


Sens.  14:1119-1131. 

Vidal,  A.  and  A.  Perrier  (1990)  Irrigation  monitoring  by  following  the  water  balance  from 
NOAA-AVHRR  thermal  infrared  data.  IEEE  Trans,  on  Geosci.  and  Rem.  Sens, 


n* 


*n 

ON 

I 

ON 

Tf 

ON 

OO 

<N 


FIGURES  AND  TABLES 


Figure  2.  Landsat  TM  image  covering  WGEW  and  the  area  immediately  surrounding.  The 
solid  lines  delimit  the  boundaries  of  the  WGEW  and  RIP  targets.  The  flight  paths  of 


VO 


Figure  5.  An  investigation  of  the  effect  of  emissivitv  and  sensor  calibration  on  aggregation 


sensible  heat  flux  (H)  using  thermal  scanner  data.  Values  of  H were  calculated  using 
Eqs.  [10]  - [12]  and  thermal  scanner  measurements  of  Tk  for  Kendall  (K)  and  Lucky 
Hills  (LH)  on  DOY  209  and  DOY  221 . was  assumed  to  be  constant  at  both  scales. 
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Figure  9.  An  investigation  of  the  effect  of  stability  corrections  of  r/  and  site-specific 


0.01  m at  MF3-MF7  and  0.04  m at  MF1,  MF2  and  MF8  for  the  bar  labeled 
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right.  "Std.  Dev."  is  one  standard  deviation  from  the  mean,  where  n=  143,1 12  in  all 
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Table  1.  General  description  of  instruments  deployed  during  the  Monsoon'90 


Table  2.  Summary  of  statistics  associated  with  histograms  of  thermal  scanner  data  (acquired  with  the  Cessna 
aircraft  flying  at  0.09  km  to  3.35  km  above  ground  level)  presented  in  Figure  4.  "Skew"  is  the  3rd 
moment  about  the  mean,  where  a larger  absolute  value  represents  a more  pronounced  skew  and  a 
negative  number  indicates  a skew  to  the  right.  "Std.  Dev."  is  one  standard  deviation  from  the  mean. 
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Figure  3.  An  illustration  of  the  variability  of  the  yoke-based  (0.3  m resolution)  and  TM  spectral  data  (120  m 


Figure  4.  Comparisons  of  histograms  of  thermal  scanner  data  to  illustrate  the  variability  of  Tr-T#  values  due 
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Figure  5.  An  investigation  of  the  effect  of  emissivitv  and  sensor  calibration  on  aggregation  of  surface 


ps6ej0AV  ueiji  p0jB|no|BQ 


I ■O 
Jo  £ 
no  tg 

D 

o 

to  (0 

* o 

c 

<D 

o sz 
- * «- 
TJ 
Q) 
D) 

jin  co 


O ' 
CO 

1- 


vo 


Figure  6.  An  investigation  of  the  effect  of  the  r„'  stability  correction  on  aggregation  of  sensible  heat  flux  (H) 


Metflux  Site 


Figure  7.  An  investigation  of  the  effect  of  the  r J stability  correction  on  aggregation  of  sensible  heat  flux  (H) 
using  thermal  scanner  data.  Values  of  H were  calculated  using  Eqs.  [10]  - [12]  and  thermal  scanner 
measurements  of  Tk  for  Kendall  (K)  and  Lucky  Hills  (LH)  on  DOY  209  and  DOY  22 1„  z was 
assumed  to  be  constant  at  both  scales,  though  different  for  each  site.  * 

Thermal  Scanner  (0.3  m resolution) 
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Figure  8.  Comparison  of  H values  computed  for  a continuum  of  (TVTJ  values  DOY  252  at  10:30  at  Kendall 
(^=0.01  m)  and  Lucky  Hills  (v=0.04  m)  for  stable  and  unstable  conditions.  These  data  illustrate  the 
non-linearity  of  the  relation  between  (Tk-TJ  and  H,  and  the  differences  in  the  relation  due  to 
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Figure  10.  An  investigation  of  the  effect  of  stability  corrections  of  r/  and  site-specific  differences  in  zj  on 
estimates  of  sensible  heat  flux  (H)  using  NS001  data.  Values  of  H were  calculated  using  Eqs.  [10]  - 
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GCM  first  level.  Each  tile  has  its  own  prognostic  variables.  The  authors  argued  that  the 
tile  approach  (the  LEAF  and  MOSAIC  surface  schemes)  are  to  some  extent  simpler  than  the 
classical  ’mixture’  big  leaf  model.  For  instance  the  resistance  network  in  Leaf  is  simpler  than 
in  SiB.  An  other  advantage  is  that  the  tile  approach  can  provide  additional  information  on 
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cover.  Figure  2 gives  the  daily  course  of  the  sensible,  latent  and  soil  heat  fluxes.  The  sensible 
heat  flux  was  about  two  times  lower  in  the  crop  area  than  in  the  forest  (Noilhan  et  at.  1991). 
On  the  other  hand,  the  evaporation  flux  was  much  less  scattered.  The  decrease  of  the  forest 
transpiration  was  interpretated  as  the  result  of  stomata!  control  due  to  increasing  atmospheric 
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planetary  boundary  layer  is  reduced  by  the  advection  of  moist  air.  On  the  other  hand,  the 
latent  heat  flux  over  the  western  Tomelloso  area  is  increasing  with  height,  due  to  the  strong 
entrainment  of  dry  air  at  the  top  of  the  boundary  layer.  The  strong  value  of  entramement  is 
associated  with  the  rapidly  growing  planetary  boundary  layer  and  the  large  discontinuity  of 
specific  humidity  at  the  top  of  the  boundary  layer.  Such  quite  complicated  processes  can  be 


clearly  exhibited  only  through  numerical  mesoscale  modelling  since  no  experimental  program 
will  ever  be  able  to  gather  such  a wealth  of  information. 
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calibration  tests  with  ISBA.  The  bi-modal  shape  of  the  pdf  of  surface  conductance  for  the  Hapex 
area  characterizes  the  partitioning  of  the  area  into  the  pine  forest  and  the  agricultural  zone.  On 


(not  shown),  with  the  exception  of  the  bare  soil  evaporation  component  in  the  case  of  large 
variability  of  soil  texture  (e.g.,  juxtaposition  of  sandy  and  loamy  soils  in  the  Hapex  experimental 
area).  Indeed,  the  effects  of  non-linearity  seem  more  pronounced  in  the  cases  of  soil  properties 
than  for  vegetation. 


2.2  Subgrid  precipitation  and  interception  loss 

Subgrid  precipitations  and  its  relation  with  interception  loss  is  an  important  issue  in  GCM 
modelling.  It  is  well  known  that  the  partitioning  of  evaporation  between  wet  and  dry  canopy 
is  incorrect  if  the  precipitation  is  uniformly  spread  in  the  grid  box.  Parameterizations  of  the 
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potential  rate.  If  the  sub-grid  distribution  of  rain  is  not  taken  into  account  in  the  ID  model, 
the  evaporation  flux  is  clearly  overestimated  because  the  interception  reservoir  is  maintained  at 
its  saturated  value.  The  reduction  of  the  fractional  wetted  area  has  the  effect  of  bringing  the 
rate  of  evaporation  down  to  levels  given  by  the  3D  model.  Additionally,  Figure  9 shows  that  the 
’GCM’  effective  evaporation  of  intercepted  rain  is  significantly  improved  if  the  daily  cycle  of  fj. 


is  applied  in  the  aggregation  scheme,  instead  of  using  a constant  value  as  usually  done  in  many 
GCMs 
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Figure  8:  Diurnal  variation  of  the  fractional  area  n of  a large  mesoscale  sub-domain  covered 
by  precipitation.  The  fractional  area  p is  computed  from  the  3D  model. 

Figure  9:  Diurnal  variation  of  modelled  evaporation  flux  of  intercepted  rainfall  on  5 June 
1986:  (solid  circles  and  vertical  bars)  3D  model;  (long-dashed  line)  ID  test  without  subgrid 


variability  of  interception;  (solid  line  ) ID  test  including  the  subgrid  distribution  of  rainfall  fj. ; 
(short  dashed-line)  ID  test  with  /z  constant  during  the  day. 
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During  the  last  decade  or  so,  procedures  to  represent  landscape  heterogeneity  within  a, 
larger- scale  area,  such  as  a general  circulation  model  (GCM)  grid  cell,  have  focused  on  a 
summation  of  surface  fluxes  within  that  area,  proportionally  weighted  by  the  fluxes  from 
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et  al.  1992). 

‘We  define  mesoscale  using  the  definition  presented  in  Pielke  (1984).  Mesoscale  systems  can  be  accurately 
described  by  the  hydrostatic  equations  of  motion,  while  the  turbulent  motions  are  nonhydrostatic  features. 

features  differ  from  larger  scale  systems  in  that  the  horizontal  winds  are  substantially  out  of 
gradient  wind  balance  even  above  the  planetary  boundary  layer. 


Therefore,  it  should  be  possible  to  parameterize  the  influence  of  the  mesoscale  fluxes 
on  larger  scale  atmospheric  structures.  This  is  the  purpose  of  this  paper.  Claussen  (1991), 
for  example,  has  discussed  how  land  surface  variability  influences  spatially  averaged'  fluxes 
and  spatially  averaged  vertical  gradients,  although  his  analysis  is  limited  to  the  situation 
where  the  boundary  layer  becomes  homogeneous  above  a “blending  height  (Wierenga, 


mogenize  the  atmosphere  a short  distance  above  the  surface  (such  that  horizontal  gradients 
in  boundary  layer  structure  do  not  develop)  when: 


Kx  TFi(tm)  = Order  (1)  (1) 

where  R0  is  the  Rossby  radius  of  deformation,  T is  proportional  to  the  inertial  period  and 
a period  related  to  deceleration  of  horizontal  flow,  m = 2Rq/Lx,  and  F\(tm)  is  the  time 


to  parameterize  them)  become  insignificant  when  the  horizontal  scale  of  the  flow  is  on  the 
same  order  as  the  vertical  scale  (i.e.,  R0/m  S ho).  Avissar  and  Pielke  (1989)  assumed 
this  result,  but  without  theoretical  justification,  in  developing  a mosaic  one-dimensional 
representation  of  land  surface  effects. 


Using  mixed  layer  scaling  (e.g.,  see  Garratt,  1992),  the  horizontal  exchange  coefficient 
can  be  estimated  as: 


| 

ai 

**? 


with  landscape  patchiness,  we  first  have  explored  the  predictability  of  these  flows.  The 
predictability  measure  that  we  apply  is  the  signal-over-noise  ratio,  r.  The  signal  is  defined 
to  be  the  domain-averaged  root  mean  square  (rms)  difference  of  a variable  at  time  t and  at 
the  initial  time,  and  the  noise  is  defined  to  be  domain-averaged  rms  difference  of  a variable 


in  the  control  simulation  and  in  the  perturbed  simulation  at  time  t.  The  domain  average  is 
obtained  from  a vertical  and  horizontal  summation  of  each  grid  point  in  the  domain.  In  order 
to  estimate  the  overall  predictability  during  the  day,  we  have  also  computed  the  mean  value 
r,  of  the  signal-over-noise  ratio  r(t)  from  t — 3 to  11  hours.  When  f(a)  > 1 the  variable  a is 
defined  to  be  predictable  during  the  day,  otherwise  the  variable  a is  unpredictable.  Figures 
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We  expect  that  mesoscale  flaxes  would  vary  in  a nonlinear  fashion  when  these  parameters 
vary,  but  would  be  a well-defined  function  of  these  effects.  For  example,  mesoscale  effects 
would  be  expected  to  dominate  for  intermediate  values  of  Lx  (Xian  and  Pielke,  1991;  Avissar 


and  Chen,  1993).  For  very  small  values  of  Ls,  the  atmosphere  is  quickly  homogenized 
immediately  above  the  surface  such  that  mesoscale  circulations  do  not  develop  (e.g.,  Hadfield 
et  al.,  1991;  1992).  For  very  large  values  of  Lx,  mesoscale  circulations  would  develop  only 
near  the  patch  edge  with  the  interiors  dominated  by  turbulent  fluxes.  This  paper  places 
quantitative  limits  on  the  influence  of  these  parameters  on  the  flux. 


to  model  grids  with  A*  = 2 km,  and  small-scale  (or  turbulent)  fluxes , which  are  defined  as 
the  subgrid  parameterized  fluxes  in  this  model  grid,  the  general  formula  is 


jo!  + 03A1  + a3Af  + a4 A4  + a5A|  + a*\e  + a7\%  + a8  (£)  + 09  (%) 
[&i  + 62 Ai  + 63A?  + 64A4  + M*  + Me  + M2  + &s  (fir)  + 69  (ft)  ] t 


,re  fa  is  a quadratic  polynomial  as  expressed  in  Eq.  (6).  The  values  of  the  coefficients 
b,-,  and  c,-  in  Eqs.  (6)  and  (7)  are  given  in  Table  6. 

For  the  total  momentum  flux  scaled  by  its  surface  value  F , the  formula  is 


near  the  surface  the  turbulent  fluxes  are  dominant;  however,  in  the  upper  portion  of  the 
boundary  layer  the  mesoscale  fluxes  (especially  those  of  sensible  heat  and  momentum)  are 
dominant.  This  conclusion  is  consistent  with  the  interpretation  of  Smith  et  al.  (1994)  and 
others  who  participated  in  the  FIFE  study,  in  that  mesoscale  fluxes  did  occur.  This  is  also 


consistent  with  the  estimate  of  FIFE-area  surface  fluxes  obtained  from  FIFE-area  averaged 
variables  which  are  essentially  identical  to  the  linear  averages  of  individual  point  specific 
fluxes,  as  reported  in  Sellers  et  al.  (1992b),  because  the  total  fluxes  are  always  dominated 
by  turbulence,  and  the  mesoscale  ascent  is  small  near  the  surface. 
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Table  9:  Same  as  Table  8 except  X*  = 20  km. 
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Table  2:  The  coefficients  a,-  used  in  Equation  (2)  for  mesoscale  [denoted  by  superscript 
“ ' ”]  and  small-scale  (or  turbulent)  [denoted  by  superscript  “ " ”]  sensible  heat,  moisture, 
and  momentum  fluxes  scaled  by  their  corresponding  total  (or  turbulent)  surface  values. 
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Table  3:  Tie  coefficients  6,-  used  in  Equation  (2)  for  mesoscale  [denoted  by  superscript 
“ ' ”]  and  small-scale  (or  turbulent)  [denoted  by  superscript  “ " ”]  sensible  heat,  moisture, 
and  momentum  fluxes  scaled  by  their  corresponding  total  (or  turbulent)  surface  values. 
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Table  4:  The  coefficients  a used  in  Equations  (2)  and  (3)  for  mesoscale  [denoted  by 
superscript  “ ' ”]  and  small-scale  (or  turbulent)  [denoted  by  superscript  “ " ”]  sensible  heat, 
moisture,  and  momentum  fluxes  scaled  by  their  corresponding  total  (or  turbulent)  surface 
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Table  5:  The  coefficients  o,-,  6,,  and  c,-  nsed  in  Equations  (4)  and  (5)  for  total  sensible 
heat  flux  scaled  by  its  surface  value. 
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Table  6:  The  coefficients  a,,  and  c,-  used  in  Equations  (5)  and  (6)  for  total  moisture 
flux  scaled  by  its  surface  value. 
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Table  7:  Tbe  values  of  the  coefficients  a„  c,-,  and  d,-  used  in  Equations  (7)  and  (8) 
for  total  momentum  flux  scaled  by  its  surface  value. 
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Table  8:  Tbe  ratio  of  mesoscale  to  turbulent  fluxes  of  heat,  moisture,  and  momentum 
as  a function  of  dimensionless  height,  z/z,-,  for  the  FIFE  study  area.  The  parameter  Lx  = 
10  km  and  other  parameters  are  given  in  the  text. 
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Table  9:  Same  as  Table  8 except  Ls  = 20  km. 
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The  Influence  of  Topography  on  Meteorological  Variables  and  Surface-Atmosphere 
Interactions 


scalars  (especially  sensible  and  latent  heat)  over  low  to  moderate  hills,  both  without  and  with 
flow  separation.  The  treatment  of  the  wind  field  is  based  on  a physical  division  of  the  flow  into 
inner,  outer  and  wake  regions.  This  framework  is  used  to  discuss  modelling  approaches, 
measurements  of  both  the  mean  flow  and  the  turbulence,  special  processes  in  the  wake  region, 
and  the  effects  of  stratification. 


Scalar  transfer  and  the  surface  energy  balance  are  discussed  first  in  terms  of  a linear 
theory  which  accommodates  radiative,  aerodynamic,  and  elevation  perturbations  on  the  surface 
energy  balance  over  hills.  The  linear  theory  provides  general  scales  for  determining  the  relative 
magnitudes  of  these  perturbations,  but  its  quantitative  predictions  are  limited  to  low  slopes 
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Common  observation  shows  that  topography  has  a profound  influence  on  the  climate  near 
the  ground,  and  thereby  on  exchanges  of  energy,  water  and  momentum  between  the  land  surface 
and  the  atmosphere.  Hill  crests  are  windier  than  valleys;  hillslopes  facing  the  equator  are  warmer 
than  dnncc  farina  the  noles:  toooeraDhv  strongly  modulates  cloudiness  and  precipitation;  it  is 
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(7)  Land-atmosphere  exchanges:  All  the  above  processes  influence  the  land-atmosphere 
fluxes  of  energy,  water,  carbon  and  momentum.  For  example,  let  us  consider  the  case  of 
evaporation  (a  function  of  radiation,  saturation  deficit,  wind  speed  and  surface  resistance)  over 
terrain  covered  by  long,  parallel,  sinusoidal  ridges.  For  simplicity,  we  imagine  that  the  ridgelines 
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energy  balance  (SEB)  in  both  local  and  spatially  averaged 
not  exclusively,  on  the  case  of  low  hills  as  defined  above. 


2.  WIND  AND  TURBULENCE  OVER  HILLS 
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length  L,  as  sketched  in  Figure  3a.  The  dP/dx  term  in  Equation  (1)  therefore  perturbs  the 


This  is  a layer  much  shallower  than  L,  because  uJUb  is  small  (typically  0.05  to  0.1).  The  outer 
region  is  the  part  of  the  flow  where  only  the  pressure  gradient  is  important  in  modifying  the 
velocity  field,  while  the  inner  region  is  the  shallower  layer  where  the  velocity  field  is  modified 
by  both  pressure  and  stress  gradients. 


The  depth  / of  the  inner  region  is  usually  regarded  as  a global  parameter  for  flow  over  a 
hill,  although  the  inner  region  as  defined  above  is  not  of  uniform  depth  (Figure  3b).  Jackson  and 
Hunt  (1975)  estimated  / as  the  height  at  which  the  perturbation  stress  divergence  and  advection 
terms  in  Equation  (1)  are  of  the  same  order  of  magnitude.  For  neutral  flow  with  a logarithmic 
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independent  of  changes  in  the  turbulence  because  the  effect  on  the  mean  flow  of  changes  in 
turbulence  is  confined  to  the  inner  region.  In  the  outer  region,  the  equations  of  motion  and 
continuity  can  be  cast  in  perturbation  form,  written  in  the  inviscid  approximation  (by  neglecting 


terms  involving  turbulent  stresses),  and  linearised  (by  neglecting  perturbation  product  terms). 
The  flow  is  then  governed  by  a single  equation  for  W,  which  in  the  two-dimensional  case  is  . 
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neutral  surface  layer,  Te  is  of  order  z/u,.  In  the  outer  region,  Te  exceeds  the  advection  time  over 


acceleration,  shear  and  curvature.  This  response  can  be  determined  from  linear  rapid-distortion 
theory"  (Hunt  and  Carruthers  1990),  in  which  the  effects  of  mean  flow  changes  on  the  turbulence 
are  calculated  from  linear  equations,  ignoring  the  nonlinear  turbulence-turbulence  interactions 


because  they  require  time  scales  > Te  to  take  effect. 
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hill;  (c)  linear  equations  are  formed  for  the  perturbation  components,  neglecting  perturbation 
product  terms  relative  to  linear  terms,  which  is  a good  approximation  if  the  perturbations  are 
small  enough;  (d)  the  linearised  equations  are  then  solved  in  the  various  layers,  with 
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incorporated  directly  into  the  coordinates,  with  the  mean  vertical  wind  (W)  being  replaced  as  a 
flow  variable  by  streamline  radius  of  curvature  (/?).  The  dynamical  effects  of  streamline 
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Figure  4 summarises  some  of  the  main  features  of  the  mean  velocity  field  over  a two- 
dimensional  (2D)  ridge  normal  to  the  mean  wind,  and  over  an  axisymmetric  three-dimensional 
(3D)  hill.  At  the  upwind  foot  of  the  2D  ridge  (but  not  the  3D  hill),  a slight  deceleration  is 
observed  in  the  near-surface  (inner-region)  flow.  This  is  in  response  to  the  decelerating  (positive) 


are  the  height  zmnT  of  the  maximum  in 


maximum.  The  linear  theory  of  HLR  predicts  that  * U3,  where  / is  the  inner-region  depth. 


where  t| (x,z,z0)  is  a function  of  order  1 , dependent  on  the  precise  hill  shape  and  the  upwind 


roughness  length  Zq  (details  are  given  in  HLR).  As  a test  of  the  HLR  prediction.  Figure  6 shows 
a transect  with  x of  measured  speedups  at  height  2 m on  Cooper’s  Ridge,  New  South  Wales 
(Coppin  et  al.,  1994).  The  HLR  prediction  is  shown  both  for  the  actual  hill  shape  and  also  for  the 
"Witch  of  Agnesi"  shape,  a commonly  used  analytical  shape  for  a 2D  ridge,  given  by 


where  Z(x)  is  the  local  elevation.  This  hill  shape  has  convenient  mathematical  properties,  and  is 
the  shape  sketched  in  Figures  2 and  4.  The  predictions  in  Figure  6 for  the  two  shapes  are  fairly 
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normalised  with  Z.  Figure  7,  like  Figure  5,  is  a synthesis  of  findings  from  the  field  experiments 
in  Table  1,  together  with  wind  tunnel  data.  The  rather  complex  behaviour  of  the  profiles  shown 
in  Figure  7 occurs  as  the  turbulence  strikes  a balance  between  the  competing  demands  of 
equilibrating  with  changes  in  the  mean  flow  at  z « / (where  the  turbulence  time  scale  Te  « 
UUb,  the  travel  time  over  the  hill)  and  responding  "rapidly"  to  shear,  acceleration  and  curvature 


inner  region.  The  near-surface  region  of  local  equilibrium  is  now  replaced  by  a plant  canopy 
flow,  modified  from  its  usual,  level-terrain  form  (e.g.  Raupach,  1988)  by  flow  acceleration  and 


curvature.  Finnigan  and  Brunet  (1994)  made  a wind  tunnel  study  of  a model  canopy  (height  hc ) 
covering  a 2D  ridge  with  height  H = 3 hc  and  length  L = 8.4/ic,  finding  a pronounced  variation 


side,  0 , exceeds  the  values  given  in 


Whether  a separation  region  is  formed  or  not,  a wake  region  of  reduced  mean  velocity 
and  enhanced  turbulence  extends  for  a considerable  distance  downwind  of  the  hill.  The  properties 


atmosphere,  defined  by  N = [ (g/0) (de/8z)]in  , with  g the  gravitational  acceleration  and 


mean  potential  temperature;  N is  the  angular  frequency  of  vertical  oscillation  of  an  air  parcel 
under  the  restoring  force  of  buoyancy.  The  Froude  numbers  FL  and  FH,  based  respectively  on  hill 

length  and  height,  are  infinite  in  a neutral  ambient  atmosphere  (N  - 0),  and  decrease  towards 
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substantially  revised  physical  view  of  this  concept  was  presented  by  Smith  (1990),  but  the  simple 
formula  compares  surprisingly  well  with  towing  tank  data  (Snyder  et  al.,  1985). 
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HLR  for  the  velocity  field.  However,  a key  distinction  between  the  scalar  and  velocity  fields  is 
evident  by  comparing  the  momentum  conservation  equation,  (1),  with  the  corresponding  scalar 


the  perturbations  in  U,  W and  K,  they  obtained  and  solved  a Unear  equation  for  perturbations  in 
5.  The  solution  shows  that  5 perturbations  arise  through  three  aerodynamic  mechanisms:  (a) 

convergence  and  divergence  of  mean  streamlines;  (b)  vertical  stress  gradients  dx/dz  , which 


because  changes  in  the  9 field  affect  the  surface  flux  of  Q (via  the  saturation  deficit),  and  vice 
versa.  Before  the  problem  can  be  analysed  for  each  scalar  individually,  these  boundary 


conditions  must  be  decoupled.  Following  McNaughton  (1976),  a general  way  to  do  this  is  to 
define  two  new  scalars,  each  a linear  combination  of  0 and  Q: 
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into  thermal  storage,  and  the  subscript  1 denoting  a background  value);  and  a dimensionless 
elevation  parameter  Pelev,  proportional  to  the  adiabatic  change  in  saturation  deficit  through  height 


vaporative  fraction  for  equilibrium  evaporation,  about  0.68  at  20  C (Slatyer  and  McDroy,  1961; 
rhom,  1975;  Raupach,  1994).  These  scaling  factors  imply  that  radiative  and  aerodynamic  effects 
ire  both  proportional  to  hill  slope  H/L,  whereas  elevation  effects  are  proportional  to  hill  height 
H.  Also,  aerodynamic  effects  make  a positive  contribution  to  the  latent  heat  flux  perturbation 


c8 


J§ 


g rr 

3 u 


8.  C 


% & 

•9  £ 

<U 

I M 

W)  *cd 

• pH  f-”H 

‘53 

& s 

0 

1 9 

5C  h* 

g * 

£ 5 


06  « •* 


SJ  cd 


c«  X) 

G * 


a kT 

a ^ 


§ 

£ 

« 

w 

TO 

< 

(Q 

o 

CO 

■i-H 

X5 

p 

< 

12 

c3 

G 

3 

O 

Vi 

cd 

P 

u 

P 

r/\ 

O 

& 

<8 

c3 

o 

<N 

g 

<o 

p 

p 

8 

m 

G 

00 

VO 

W) 

•s 

o 

H 

It 

ll 

A 

Vi 

r-H 

•a 

C 

1 

5 

e 

p 

T) 

Vi 

§ 

1 

'w' 

Si 

iH 

3 

Vi 

3 

pS 

a 

* 

o 

13 

> 

c W 
a U 

* £ 


VO  D 
w U 

S " 
I & 

3 5* 

CT1  ”5 

B 3 

p-H 

3 -s 


<D  CO 

§ M 


■9  H. 

1 1 
•g  s 

G T3 


^ \G 

3 « 


u O 

I 3 


.&  * 


£?  -o 

’S 

£ .S 


& % 

o 2 


X 

3 

B 

53 

31 

•w 

!3 

Cd 

pG 

G 

P 

a 

G 

O 

C/3 

G 

C/3 

c 

P 

C/5 

T3 

o 

•9 

■8  ^ 


<D  cfl  5 

•a  b .g 

e u g 

<D  ^ 2 

IS  I 

to  o 

£ q*  +-» 


O cd 

p a 


p P 

*3  .2 

a i 


« < 
I 3 s 


2 

h? 

w 

s 

S 

c/3 

-*-» 

O 

P 

43 

4— • 

G 

♦ 

cd 

O 

53 

♦ *H 

G 

< 

u 

P 

4H 

OX) 

• pH 

T3 

*3 

C/3 

2 

C/3 

9 

1 1 

<d  <u 

p a 


p 


S3  < Z 

S jS  5 

0 3 m 

r£  Cd  p 

to  w > 

« C O, 

D O 

> P P 

1 a sp 


*2  .2 
S -o 

W i) 


"Joe 

S.  ■s  I 

2 8 c 

u 1 I 

S £ I 


u k 

9 g 

f \ H 


td  -fl 

> .pH 


>*  Vi 
Vi  o 


' W)  O 

«*  G G 


> 3 


5 t> 

s 8. 

<D  on 


S3  «a 

O *=i 


0)  P 

Q<  C 


perturbations  in  available  energy  is  exactly  zero,  while  the  spatial  averages  of  the  FE  and  F, 
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Nonlinear  Effects:  an  Experiment  on  Scalar  Transfer  over  Hills 


The  consequences  of  separation  and  other  nonlinear  processes  in  the  wake  region  are 
difficult  to  assess  theoretically,  so  we  must  make  best  use  of  experimental  techniques.  Wind 
tunnel  measurements  have  a great  deal  to  offer  in  studying  these  complex  phenomena,  because 
the  spatial  density  of  measurements  required  is  very  large  and  because  field  experiments  are 
subject  to  natural  vagaries  in  the  weather.  Of  course,  field  experiments  are  crucial  and  have  a 
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Several  checks  were  first  carried  out  to  establish  the  momentum  and  heat  transfer 
properties  of  the  three  roughnesses  on  a horizontal  surface  in  the  absence  of  a hill.  These  showed 
that  the  measured  roughness  lengths  for  momentum  and  heat  for  all  surfaces  were  consistent  with 
available  published  data. 


To  examine  the  perturbations  in  heat  transfer  caused  by  a hill,  we  use  the  Stanton  number 
heat  transfer  coefficient  from  the  surface  to  the  free  stream,  CH(x): 


distance  x reflected  changes  in  the  surface  temperature  90  caused  by  variations  in  aerodynamic 


downwind  slope.  More  detailed  scrutiny  indicates  several  points,  all  of  which  were 
experimentally  reproducible.  Treating  first  the  steep  hill  2 (where  separation  always  occurred), 

the  region  of  increased  C„l  in  the  separated  wake  region  extended  about  2 L downwind  from 
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For  flat  terrain,  the  influence  of  heterogeneity  on  the  regional  SEB  has  been  the  subject 
of  much  recent  work.  Raupach  and  Finmgan  (1994,  henceforth  RF)  tested  a scale-insensitivity 


hypothesis,  that  the  spatially  averaged  SEB  of  a patchwork-quilt  surface  is  insensitive  to  the  scale 
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Here,  we  apply  a similar  approach  to  infer  the  effect  of  topography  on  the  regional  SEB. 
The  approach  is  to  scale  up  the  above  measurements  of  topographically  induced  heat  transfer 
perturbations,  from  hills  in  a small  wind  tunnel  to  hills  1000  times  larger  in  a CBL.  The 
procedure  is: 


( 1 ) We  assume  a vegetated  land  surface  with  the  simple  properties  specified  in  Table  4.  This 
surface  is  either  flat  or  is  "crinkled"  with  a series  of  ridges,  lateral  to  the  wind,  with  the  shape  of 
hill  2 in  the  wind  tunnel  experiment  described  above,  but  1000  times  larger;  that  is,  L = 125  m 
and  H = 43  m.  The  ridges  are  8 L (1000  m)  apart,  a separation  sufficient  to  ignore  interactions 


speed  in  the  wind  tunnel  experiment  (Figure  12).  This  spatially  dependent  resistance  includes  all 


the  aerodynamic  effects  occurring  around  the  hill,  as  it  is  obtained  experimentally  from  a 
dynamically  similar  (though  geometrically  much  smaller)  situation. 


Figure  13  shows  the  x variation  of  the  energy  fluxes,  at  0900  and  1500,  with  the  bulk 
surface  resistance  spatially  uniform  at  100  s m1.  There  are  significant  perturbations  in  FE,  FH  and 
the  available  energy  FA  = FE  + FH,  of  opposite  sign  in  the  morning  and  afternoon  as  the  sun 
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mesoscale  circulations  are  both  ignored.  If  further  work  shows  that  topographic  perturbations  on 
spatially  averaged  scalar  fluxes  are  small  even  in  the  presence  of  such  complicating  (and  highly 


subject  is  complex  and  very  wide-ranging,  and  present  efforts  are  only  the  beginning. 
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Table  1 : Some  major  field  experiments  on  flow  over  low  hills,  showing  hill  height,  H\  length, 
L;  roughness  length,  Zo;  inner  layer  depth,  /,  from  Equation  (2);  and  observed  height  of*  maximum 
speedup, 
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Table  2:  Conditions  for  separation  in  the  wake  region,  for  nonstratified  hill  flow. 
Separation  is  likely  when  maximum  hill  slope  on  the  lee  side,  0 , exceeds 


Table  4:  Surface  properties  assumed  for  calculations  of  SEB  over  repeated  ridges 
using  CBL  slab  model  and  wind  tunnel  specification  of  ra{x) 
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Figure  9:  Predictions  of  RWCH  theory  for  perturbations  in  surface  latent  heat  flux  FE,  sensible 
heat  flux  FH  and  available  energy  flux  FA.  See  text  and  RWCH  for  details. 


Figure  10:  Preliminary  wind  tunnel  experiment  on  heat  transfer  over  2D  ridges  of  varying 
steepness  and  roughness:  (a)  tunnel  arrangement;  (b)  construction  of  heated  rough  surfaces;  (o) 
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Figure  14:  Predictions  of  CBL  slab  model  with  wind  tunnel  ra(x)  for  diurnal  course  of  (a)  (F£) 


and  (b)  , for  flat  and  ridged  landscapes  with  rs  = 10,  100  and  1000  s m*1. 
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Fi2ure  3:  l a)  Schematic  diagrams  of  typical  perturbations  in  ta)  pressure  P< x 
stress  xix.z),  for  nonstratified  flow  over  a low  hill. 


Maximum  speed-up 
over  the  crest 


mean  wind,  and  a 3D  axisymmetric  hill. 


flow  over  a low  hill. 


(ui)  wBjSH 


o o o 

o o o 

CO  CD  ""t 


’J) 

D 

— i 

1 

d 

r 

X) 

a 

u> 

»cj 

p 

o 

d 

in 

X) 

* 

§ 

<u 

V 

Pi 

d 

i> 

b/> 

Jp 

'70 

TJ 

k-H 

2 

• .-4 

X2 

'✓> 

13 

Vi 

p 

V 

cx 

o 

u 

d 

o 

l) 

U 

Jp 

i— j 

*— i 

_c 

d 

M 

o 

p. 

X> 

P 

l~i 

XJ 

o 

0) 

'Vi 

1) 

CX 

,»*s 

'JO 

»n 

v_. 

X) 

K 

c 

o 

P 

• *-4 
4-1 

d 

d 

Ji 

p 

cr 

6 

W 

3 TJ 
bf)  (U 


Infrared  thermometer 


(IHM ou‘x)  tHo I (HN'x)  tH0 


J 5 


i-J 

T3  W 

u U 


2 c 

3 o 


^(srrf1)  =1000 


The  impact  of  using  area-averaged  land  surface  properties  - topography. 
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The  relationships  linking  surface  and  root-zone  soil  wetness  to  the 
soil  surface  and  canopy  transpiration  rates  are  non-linear. 


However,  simulation  results  and  observations  indicate  that  soil 
moisture  variability  decreases  significantly  as  an  area  dries  out. 


length  scales,  to  within  an  acceptable  accuracy  for  climate  mo> 


CO 


Introduction 


versions  of  their  soil-  vegetation-atmosphere-transfer  (SVAT)  models  at  GCM 
grid  scales  using  assumed  area-averaged  lower  boundary  conditions  (vegetation 
properties,  topography,  soil  moisture)  and  atmospheric  and  radiative  forcings. 
The  use  of  increasingly  sophisticated  SVATs  in  GCMs,  see  for  example  Dickinson 


(1984),  Sellers  et  al.  (1986;  1995a,  b),  Bonan  (1994),  has  almost  always  resulted  in 
improved  model  performance  and  the  simulation  of  more  realistic  continental 
climatologies,  see  for  example  Sato  et  al.  (1989),  Randall  et  al.  (1995),  Dickinson 
and  Henderson-Sellers  (1988),  Nobre  et  al.  (1991).  However,  in  parallel  with  this 


Upscale  integration  of  models:  Experiments  were  designed  to  test  the  soil- 
plant-atmosphere  models  developed  by  biometeorologists  for  small-scale 


from  the  First  ISLSCP  Field  Experiment  (FIFE)  to  explore  the  impact  of  area- 
averaging initial /boundary  conditions  of  topography,  vegetation  condition  and 
soil  moisture  on  the  calculate  fields  of  radiation,  heat  and  water  vapor,  see  Sellers 
et  al.  (1992a).  In  the  process,  methods  of  defining  surface  properties  from 
satellite  data  were  developed  and  tested. 


Theoretical  Consideration  in  'Upscaling'  Soil-Vegetation- Atmosphere 
Transfer  (SVAT)  Models 


SiB  was  modified  in  Sellers  et  al.  (1992b)  to  run  independently  of  a GCM  with  the 
reauired  atmospheric  forcings  provided  by  near-surface  observations  of 


tf) 


an  extension  of  the  Penman-Monteith  equation,  see  Monteith  (1973).  Like  just 
about  all  SVATs  in  use  today,  SiB  has  been  demonstrated  to  provide  realistic 
local-scale  simulations  of  the  surface  sensible  and  latent  heat  fluxes  for  many 
different  vegetation  types,  see  for  example  Sato  et  al.  (1989),  Sellers  et  al.  (1992b). 


Thus,  if  we  have  a model  for  Q as  a function  of  xi  which  we  believe  works  well  at 
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<Q>I  = estimate  of  area-averaged  value  of 
surface  flux  provided  by  numerical 


The  angle-brackets  < >'  denote  'area-average'.  Figure  2 (left  hand  side)  shows 
how  a numerical  calculation  may  be  done  to  yield  <Q>i  by  dividing  up  the 
domain  A into  many  small-scale  segments  for  each  of  which  (1)  can  be  assumed 
to  be  valid.  In  this  approach,  the  fields  of  xi  are  discretized  over  the  domain  so 


Figure  2 (right  hand  side)  shows  how  this  use  of  (3)  reduces  the  tedious 
numerical  integration  of  (2)  to  a single  calculation  but,  as  f (xj)  becomes  more 

A 

non-linear,  the  difference  between  the  'true'  value  of  <Q>i  as  given  by  (2),  and 
the  area-averaged  estimate  <Q>m  provided  by  (3)  will  increase. 
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within  the  domain  depends  on  slope  and  aspect,  see  Figure  4.  In  this  study,  w 
adapted  the  SiB  radiation  submodel  to  take  account  of  illumination  geometry. 
The  modification  is  simple  and  does  not  take  into  account  the  effects  of 
topographic  shadowing  by  other  pixel  elements  in  the  scene. 


[S4-b,v  (l"ab,v)  + S'lb/n(l'ab/n)]Pb 


0> 

£ 


xJ  u 
<U 

Vh  fV 

r*  t/i 

C (tf 

S *3 

•°  i 

4-» 

u <u 
2 Oh 


<U  M 

45 

s | 

<U  •> 

3 I 

> ^ 
o>  <u 

e e 

o o 

nd  X* 


o o 

e g 

o o 

*43  *43 
u u 
(ts  as 

I-* 

Ph  PL, 


43  <u 


s 


a ^ g u 

g>  s s s * 

is  8.8.5. 

£ o £ ,2  <£> 

c n ui  in  *•* 


6 -5 

Q <1 1 


T3  ih 
CO.  CQ. 


N > 8 


r!  tTj 

.5  0) 
cu 

4-*  4-* 

tn 

m XI 
C 

0)  R 


fourth  terms  cover  the  absorption  of  radiation  reflected  from  the  surrounding 


landscape  which  would  be  exposed  to  a sloping  pixel  from  the  horizon 
downwards.  Note  that  we  assume  that  the  landscape  has  the  same  reflectance 
properties  as  the  subject  pixel  and  that  there  is  no  shadowing  of  the  subject  pixe 
by  distant  topography.  (A  full-up  radiative  transfer  treatment  of  these  issues 
would  be  very  complex  and  computationally  expensive). 
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The  effects  of  topography  on  the  radiation  budget  are  encapsulated  in  the  factors 
(3b/  Pd  and  (3r,  see  (5b,  c,  d).  Generally  speaking,  deviations  from  a flat  surface 
will  lead  to  an  increase  in  the  absorption  and  emission  of  longwave  radiation 
with  a probably  small  net  gain  for  the  surface  as  the  fraction  of  cold  sky  viewed 
by  the  surface  decreases.  The  impact  of  topography  on  the  interception  of  diffuse 


radiation  is  fairly  small  as  the  teams  multiplied  by  pd  and  pr  tend  to  decrease 
slowly  from  the  solution  for  a level  surface  as  V increases.  The  absorption  of 
direction  beam  radiation  is  almost  linearly  related  to  V through  Tan  V m pb-^ 
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:unction  of  FPAR  and  vegetation  physiology  only.  Figure  (5a)  shows 
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is  because  the  transpiration  formulation  approximately  satisfies  the  condition 
represented  by  (4a) 
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the  leaf  water  potential  stress  factor  in  (7).  An  expression  for  f (V*)  and  the  other 
stress  factors,  f(5e)  and  f(T)  were  obtained  from  curve  fits  to  the  data  of  Verma  et 
al.  (1993)  such  that: 
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(1992b).  In  Sellers  et  al.  (1992b)  a series  of  inverse-mode  runs  with  SiB  were  used 


to  estimate  soil  evaporation  rates  and  soil  surface  resistances,  rsurf,  from  the  time- 
series  of  flux  station  data:  this  procedure  resulted  in  the  formulation  of  an 
empirical  expression  relating  rsurf  to  the  wetness  of  the  uppermost  soil  layer,  Wi 


The  FIFE-89  ’testbed'  data  set  was  selected  for  the  study  of  Sellers  et  al.  (in 
prep.)  which  focused  on  the  impacts  of  spatial  variability  on  area-averaged 
fluxes.  Much  of  the  work  reported  in  Sellers  et  al.  (in  prep.)  is  summarized  here. 


The  data  collection  effort  in  FIFE-89  was  concentrated  primarily  within  a 2 x 15 
km  strip  oriented  north-south  in  the  middle  of  the  FIFE  area.  This  so-called 
'testbed'  area  contained  most  of  the  surface  flux  stations  and  was  also  made  the 
first  priority  for  airborne  data  collection. 
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Atmospherically-corrected  SPOT  and  Landsat  TM  surface  spectral 
reflectances  at  30  x 30m,  combined  into  SR  fields. 


(vii)  Airborne  fluxes:  The  National  Research  Council  of  Canada 

contributed  the  Twin  Otter  flux  aircraft  to  both  FIFE~87  and  FIFE 


89.  The  aircraft  is  equipped  to  measure  the  turbulent  fluxes  of 
sensible  and  latent  heat,  CO2  and  momentum,  see  Desjardins  et  al. 

(1992a),  and  MacPherson  et  al.  (1992). 
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continuous  measurements  of  the  near-surface  meteorological 
variables;  temperature,  mixing  ratio,  wind  speed  (T,  q,  u); 
downwelling  shortwave,  longwave  and  net  radiation  fluxes  (Si, 
Li,  Rn);  and  precipitation  (P)  over  the  area.  Some  comparisons 
showed  that  except  for  precipitation  there  was  very  little  spatial 


The  Vp , FPAR  and  leaf  area  index  (LAI)  values  for  each  pixel  were  specified  from 
SR  fields  as  derived  from  TM  or  SPOT  data  using  equation  (8),  and  methods 
discussed  in  Sellers  et  al  (1992a,  b).  The  LAI,  canopy  greenness  fraction  and  soil 
background  reflectance  values  are  used  by  the  SiB  radiation  submodel  to 
calculate  surface  albedo  and  the  absorption  of  radiation  by  the  canopy  and 
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of  Wang  et  al  (1992).  The  soil  moisture  contents  for  the  root  zone,  W2,  and  the 
deep  soil  layer,  W3,  are  more  difficult  to  specify  directly  from  observations;  all 
we  have  are  the  neutron  probe  readings  for  the  AMS  and  flux  station  sites.  Usmg 
these  and  the  other  available  FIFE-89  data,  a study  was  made  into  which  of  the 
following  variables  are  the  best  predictors  of  W2  and  W3:  SR,  Wi  (from  the 


PBMR  data),  slope  and  elevation.  The  best  fits  were  provided  by  regression 


Generally,  model  runs  were  started  at  some  time  close  to  the  acquisition  of 
the  PBMR  data  used  to  initialize  Wi.  This  is  because  all  of  the  other  initialization 
data  are  more-or-less  continuously  available  (temperatures),  slowly  changing  in 


time  (VF/  FPAR  and  LAI)  or  constant  (soil  properties,  topography).  For  the  full 
resolution  calculations,  the  energy  and  water  balances  were  calculated  at  hourly- 
time resolution  for  each  of  the  (68  x 501  = 34068)  30  x 30m  pixels  within  the 
testbed  domain. 
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comparison  with  overlapping  2x4  km  blocks  of  footprint-adjusted  surface  flux 
estimates  produced  by  the  Twin  Otter,  see  Figure  10c  and  Desjardins  et  al. 
(1992a).  The  flux  data  of  Desjardins  et  al.  (1992a)  are  not  directly  comparable  to 
the  calculations  as  the  domains  are  different  (2x4  km  versus  2 x 2 km)  and  it  is 


known  that  estimates  of  the  surface  fluxes  provided  by  the  flux  aircraft  are 
consistently  low,  see  Desjardins  et  al.  (1992b)  and  Sellers  et  al  (1992b).  However 


the  vegetation  condition  (Vf,  albedo,  etc.)  and  surface  soil  wetness  (Wi).  The 
model  also  appears  to  reproduce  the  responses  of  the  energy  balance  to  changes 
in  soil  moisture  at  this  spatial  scale,  as  can  be  seen  by  the  range  of  EF  values  in 
Figure  10. 


Investigating  the  Effects  of  Inhomogeneitv  on  Intermediate  (2  x 15  km 
Scale  Fluxes 
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of  the  area-average  'block'  runs  where  a single  simulation  is  done  for  the  entire 
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changes  in  vegetation  condition  and  topography  are  more-or-less  linear  in  these 
variables,  so  the  impact  of  area-averaging  those  variables  is  weak.  This  result 
was  anticipated  for  the  vegetation  component  by  Sellers  et  al.  (1992b)  and  for  the 
effects  of  topography  on  the  radiation  balance  by  Dubayah  et  al.  (1990). 
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mismatch  between  the  area-integrated  (line  1)  and  area-averaged  (lines  5 and  8) 
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sensible  and  latent  heat  flux.  In  the  next  section,  we  therefore  concentrate  on  the 
influence  of  spatial  variability  in  the  soil  wetness  field  on  the  mean  surface  heat 

fluxes. 


Sensitivity  tests:  the  impact  of  soil  wetness  variability  on  the  surface  heat 
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effect  of  preserving  the  soil  moisture  variability  and  forcing  conditions  for  the 
testbed  area  while  adjusting  the  mean  soil  wetness.  (Individual  pixel  values  of  W 
were  not  permitted  to  drop  below  0.05  or  exceed  1.0  during  this  test;  the  values 
of  <W>  plotted  in  Figure  12  are  the  means  taken  over  the  W field  following  these 


bounding  adjustments  so  that  at  low  and  high  values  of  <W>,  successive 
increments  and  decrements  of  <W>  may  be  slightly  less  than  0.05).  Figure  12 
snows  that  for  the  day  216  pattern  of  soil  wetness,  the  fully-resolved  and  area- 
averaged  calculations  of  EF  are  close  to  each  other  at  high  <W>  but  diverge  as 
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decreased  in  this  test.  Are  these  patterns  conserve* 


4.2.2  Case  II:  Changes  in  soil  moisture  variability  during  a real  drvdown 


figure  13c  is  particularly  interesting:  the  spatial  variability  in  the  surface  layer 
wetness,  Gw-^  is  seen  to  decrease  as  Wi  decreases,  see  also  Figure  13e  which 

shows  the  same  data  replotted  for  clarity.  The  likely  explanation  for  this  trend  i 
that  wetter  areas  tend  to  dry  out  and/or  drain  much  quicker  than  dryer  ones 


which  has  the  effect  of  continually  narrowing  the  range  of  soil  moisture  contents 
within  the  domain  as  the  drydown  progresses.  The  results  in  Figure  13  can  be 
summarized  as  follows: 
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(iii)  As  the  domain  dries  out,  the  spatial  variability  of  the  surface  and 
root  zone  soil  wetness  decreases. 
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Figure  14  (Case  III)  compares  the  results  of  fully-resolved  and  area-averaged  SiB 
runs  starting  from  day  216  initial  conditions  but  forced  by  repeated  cycles  of  the 
day  216  meteorological  conditions  for  every  day  thereafter.  Figures  14a,  b and  c 
correspond  to  Figures  13a,  b and  c.  It  can  be  seen  that  the  time-series  of  energy 
balance  components,  EF  and  the  average  soil  wetness  calculated  by  the  two 


methods  track  each  other  very  closely  over  the  course  of  the  drydown.  This  is 
not  what  we  might  expect  from  Case  I (Figure  12)  which  indicated  that  the  two 
methods  would  produce  increasingly  divergent  surface  flux  estimates  as  W 
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case  for  days  216  through  224  shown  in  Figure  13.  The  dashed  line  shows  the 
results  of  Case  III,  the  idealized  (cyclically  repeated  day  216  meteorological 


other  closely  during  a dryd< 


Soil- Vegetation- Atmosphere  Transfer  (SVAT)  models  are  usually  developed 
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resolution  from  the  available  observations,  these  ruii-resoiunon  runs  were 
followed  by  similar  runs  in  which  one  or  more  of  the  spatially-varying  fields  was 
replaced  by  a domain-average  value.  It  was  found  that  the  use  of  area-averaged 
fields  of  topography  and  the  vegetation  index,  SR,  had  relatively  small  impacts 
on  the  mean  fluxes.  This  is  mainly  because  the  relations  linking  topography  and 
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invariant.  The  relationships  linking  the  simple  ratio  vegetation 
index  (SR),  the  canopy  conductance  parameter  (Vp)  and  the  canopy 
transpiration  flux  are  also  near-linear  and  similarly  scale-invariant 
to  first  order.  Because  of  this,  simple  area-averaging  operations  can 
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with  atmospheric  models,  i.e.  50  - 500  km?  The  results  of  this  paper  and  those  in 
Sellers  et  al  (1992b)  indicate  that  simple  area-averages  of  (moderate)  topography 
and  vegetation  parameters  are  probably  sufficient  over  a wide  range  of  spatial 
scales  provided  that  the  vegetation  type  does  not  change  radically  within  the 
domain;  the  work  of  Noilhan  et  al  (1991)  supports  this  view.  However,  the  issue 


is  not  so  straightforward  with  respect  to  soil  moisture.  To  address  it  properly,  it 
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using  area-averaged  values  of  soil  moisture,  vegetation  cover  and  topography  on 
the  calculated  fluxes  of  sensible  and  latent  heat  are  much  less  than  were 
originally  feared  by  the  authors.  This  finding  in  turn  suggests  that  inaccuracies 
associated  with  the  assumption  of  biophysical  model  scale- invariance  are  likely 


to  be  small. 
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Table  1:  Results  of  exploratory  tests  described  in  section  4.1. 
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Figure  3 Effect  of  using  area-averaged  values  of  surface  properties  xi  to 
calculate  surface  fluxes,  Q = f(xO- 
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Figure  4 Interception  of  (a)  direct  beam  radiation  and  (b)  diffuse  radiation  as 
treated  in  this  modified  version  of  SiB.  The  factor  Pb  in  equation 
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Figure  8 Schematic  showing  the  FIFE-89  testbed  area  data  sets  used  in  this 

° study.  The  figure  also  shows  the  location  of  the  testbed  area  in 


gravimetric  data  were  used  to  develop  regression 
equations  to  specify  these  soil  moisture  fields. 

(iv)  Topography:  Digital  Elevation  Model  (DEM)  from 

USGS  data,  aggregated  to  30  x 30  m resolution. 

Forcing  variables 


Near-surface  meteorology  and  downwelling  radiation 
fluxes:  Automatic  Meteorological  Stations  (AMS) 
measured  near-surface  meteorological  variables  and 
downwelling  shortwave  and  longwave  radiation 
fluxes. 


then  be  treated  as  one  big  pixel.  Any  intermediate  calculation. 


where  at  least  one  field  is  represented  at  full  resolution,  must  be 
carried  out  at  full  resolution. 
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(a)  Annual  average  Tmax  over  the  Columbia  Basin  stud} 
(b)  Annual  average  Tmin  over  study  region. 
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Abstract. 

The  effects  of  small-scale  heterogeneity  in  land  surface  characteristics  on  the  large- 
scale  fluxes  of  water  and  energy  in  the  land-atmosphere  system  has  become  k central 
focus  of  many  climatology  research  experiments.  The  acquisition  of  high  resolution  land 
surface  data  through  remote  sensing  and  intensive  land-climatology  field  experiments  (like 
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The  issue  of  ’scale  interaction’  for  land  surface-atmospheric  processes  has  emerged  as  one 
of  the  critical  unresolved  problems  for  the  parameterization  of  climate  models.  To  help 
resolve  this  issue,  the  understanding  of  the  scaling  properties  of  water  and  energy  fluxes 
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The  current  scientific  thinking  on  this  issue  is  mixed.  For  example  Sellers  et  al.  (1992) 
state  that  land- atmospheric  models  are  almost  scale  invariant,  based  on  analyses  of  data 
collected  during  the  First  ISLSCP  Field  Experiment  (FIFE);  a conclusion  also  reached 
by  Noilhan  et.  al.  (this  issue)  using  HAPEX-MOBILHY  data  and  suggested  by  Wood 
and  Lakshmi  (1993).  Counter  arguments  have  been  made  by  Avissar  and  Pielke  (1989, 
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rameters.  In  most  cases  some  terms  in  (1)  will  be  aggregated  while  others  will  be  repre- 
sented in  a distributed  manner.  The  distributed  representation  can  be  either  deterministic, 
in  which  the  actual  pattern  of  variability  is  represented  - examples  include  the  European 
Hydrological  System  model  (SHE)  (Abbott  et.  al„  1986a, b)  and  the  3-D  finite  element 
models  of  Binley  et  al.  (1989)  or  Paniconi  and  Wood  (1992);  or  statistical,  in  which  the 
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models  has  further  contributed  to  the  limited  progress. 


its  short  duration  (often  one  summer  season  or  less)  which  may  not  provide  me  range  or 
conditions  necessary  to  test  aggregation  hypotheses.  Thus,  given  the  nature  of  the  sci- 
entific debate  concerning  aggregation,  it  is  important  to  develop  a more  general  analysis. 
In  particular,  this  paper  tries  to  lay  out  a consistent  framework  for  determining  when 
land  surface  conditions  permit  aggregation  of  subgrid  variability  and  the  specification  of 


equivalent  parameters,  and  conditions  when  such  aggregation  will  fail. 

Results  presented  in  section  2.2  show  that  during  periods  of  low  atmospheric  demands 
(early  morning/late  afternoon,  Winter  periods,  etc.)  the  use  of  equivalent  parameters 
will  over  estimate  evaporation;  and  during  periods  of  high  demand  (mid-day  Summer 


necessary  for  understanding  heterogeneity  within  natural  watersheds.  Field  data  shows 
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uniform  soil  moisture  profiles.)  Thus,  given  the  atmospheric  evaporative  demand  and  the 
statistical  distribution  of  soil  moisture,  the  distribution  of  the  actual  evapotranspiration, 
Ea,  can  be  derived.  For  simple  functional  forms  the  mapping  of  z -*  d -*  Ea  can  be  done 
analytically;  and  in  any  case  it  can  be  done  through  simulation. 

Figure  1,  from  Wood  (1994),  provides  some  TOPLATS  model  simulation  results  to 


illustrate  how  the  statistical  distribution  of  the  actual  soil  evaporation  can  be  estimated 
and  by  analogy  evapotranspiration.  Results  are  provided  for  two  times  during  the  day:  an 
early  morning/late  afternoon  time  where  the  atmospheric  evaporative  demand  (potential 
rate)  is  approximately  0.35  mm/hr  and  a mid-day  time  with  an  atmospheric  evaporative 
demand  of  approximately  0.78  mm/hr.  For  surface  soil  moisture  levels  too  low  to  meet  this 
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for  surface  soil  moisture.  This  is  not  shown  here.  The  upper  left  portion  of  the  figure, 
panel  (iii),  gives  the  relationship  between  surface  soil  moisture  and  actual  evaporation  for 
the  two  times  during  the  day.  The  maximum  evaporation  rate  is  the  potential  rate,  which 
is  lower  during  the  early  morning  and  late  afternoon  (due  to  lower  net  radiation.)  For  areas 
within  the  watershed  (or  grid)  where  surface  soil  moisture  is  high,  the  actual  evaporation 


rate  is  equal  to  the  potential  rate.  For  drier  areas,  the  rate  is  lower.  The  resulting 
probability  distribution  for  the  actual  evaporation  is  shown  in  the  lower  left  portion  of 
the  figure.  This  panel  is  divided  into  two;  the  top,  panel  (iv-a),  giving  the  distribution  for 
the  time  related  to  early  morning/late  afternoon  (lower  potential  evaporation  rate)  while 
the  bottom,  panel  (iv-b),  being  related  to  mid-day  high  potential  evaporation  condition. 


d 

fn 

O 

a 


fn 

d 

03 


co 

03 


O 

£ 


03 

A 


u 

<u 

> 

03 

A 


u 

O 


a 

44 


03 

> 

03 

u 

CO 

d 

o 


o •*+ 

<h-h 

d 
o 
u 


CO 

d 

o 

HP 

d 

A 

• rH 

u 

W 


A 

rt 

4 

O 

u 

a 

-d 

v 

> 

• fH 

u 

03 

4 

03 

A 


d 

.2 

HP 

u 

03 

a 

CO 

a 


03 

fn 

0 


0 
6 

1 


& 

03 

X! 

h-> 

4 

§ 

a) 

b0 

d 

u 

* 

O 

M 

fn 

«3 

d 

d 

CO 

aj 

4l 

d 

o 

• rH 

HP 

o3 

fH 

O 

a 

> 

a; 


a 

o 


"d 

9 

p 

A 

• rH 
fH 

HP 

CO 

9* 

• rH 

*d 

>> 

1 

HP 

"d 

• r-H 

• H 

a 

"d 

1 

43 

s 

O 

fH 

O, 

d 

03 

o 

44 

o 

HP 

d 

CO 

d 

CO 

03 

»h 


A 

H 


03 

fH 

d 

o 

. rH 

u 

O 

a 

> 

03 

03 

b0 

o3 

aJ 

& 

03 

44 


B 


CO 

03 


T3 

03 

co 

d 

03 

A 

s 

u 


<u 

$ 

a) 

'bo 

{4 


03 

bO 

§ 

M 

ii 


a 

03 

d 


(4 

03 

co 

CO 

0) 


a 

.2 

HP 

u 

PS 

£ 

a 

o 

• rH 

u 

O 

& 

> 

0) 


a; 

fn 

d 


o 

6 


a 

o 


4 

a 

o 

u 

>> 

fn 

4 


a 

v 

4 

03 

1 


*h 

03 

-p 

(4 

£ 

03 

bO 

fH 

03 


a a 

.2  .S 

HP  44 

u 

d * 
.41  03 


4S 


a; 

s 


„ to 

a B 

.2  a 


o 

CO 

03 

44 

4-3 

<D 

CO 

0 

d 

U 

03 

43 

CO 

a; 

CO 

« rH 

fn 

«3 


4 

d 

O 

U 

fH 

<L> 

hp 

t! 

£ 

0) 

A 

43 

fH 

£ 


co 

d 

O 


4 

14 

O 

o 

>> 

fH 

4 

<u 

A 

H-3 

bO 

d 

• rH 
HP 

d 

03 

CO 

0) 

fH 

PU 

03 

fH 

CO 

1? 

fH 

d 


(4 

O 

u 

4 

0) 

44 

co 

fH 

03 

H-3 

ctf 

£ 

03 

A 

HP 

*4 

• *-H 

44 


O 

6 


o 

co 


s 

HP 

o 

a 


CO 

<D 

HP 

o3 

fH 

d 

o 

* 

o3 

fH 

o 

a 

> 

a; 

bO 

d 

• H 

£ 

44 


2 

3 

I 

co 

d 

M 

d 

o 

• rH 

Id 

fH 

o 

8* 

> 

a> 


d 

o 

• rH 

d 

4^ 


co 


O 

B 


o 

CO 


<D 

bO 

8 

fH 

<L> 

44 

HP 

'd 

d 

cd 

fH 

cd 

v 

d 


d 

o 

d 


a? 

& 

"d 

d 

«d 

CO 

V 

Id 

fH 

d 

o 

• rH 

"cd 

O 

& 

£ 

r. 

r— H 

bO 

d 


% 

A^ 

>> 

M 

"d 


a; 

> 

v 

u 

cd 


03 

A 

HP 

d 

03 

41 

03 

*H 

d 

HP 

03 

A 

HP 

bO 

d 

CO 

HP 

£ 

CO 

d 

d 

cd 

A 

HP 

CO 

fH 

I 

03 

Si 

"d 

d 

0) 

41 

HP 

§ 

A 


fH 

03 


Se 

t 

03 

03 

fH 

8 

CO 

d 

o 


co 

03 

bO 

bO 

d 

co 


co 

03 


03 

co 

03 

A 


co 

i! 

H 


03 

HP 

cd 

fH 

d 

o 

• rH 

Id 

fn 

o 

& 

03 

T3 


o 

fH 

fH 

03 


co 

03 

bO 


03 

A 


03 

fH 

*2 

03 

fH 

03 

A 


"d 

d 

<d 

d 

o 


g -3 

1 B 


03 

A 


"d 

d 

o 

u 

03 

44 


cd 


"d 

03 

HP 

cd 

fH 

HP 

d 

0) 

<j 

d 

o 

u 

03 

43 


O 

co 


03 

bO 

S 

fH 

03 

"d 

• rH 

cd 

A 


co 

cd 

03 

Id 


rH  0 

Id 

o | 

a 2 

f 8 

« O 

d /wN 

O fH 

03 

03  hp 

fn  03 

d d 


O 

CO 


U 

,03 

03 

fH 

cd 

0) 

d 


d 

o 

d 


cd 


CO 

d 

o 


*d 

d 

o 

u 

bO 

d 

• rH 
>> 
fH 

"d 

fH 

o 


o 

co 


d 

o 

. rH 

Id 

■ rH 

§ 

O 

HP 

03 

d 

-d 


, _ 

CO 

d 

.2 

HP 

"d 

<5 

HP 

cc3 

0^ 

CO 

03 

HP 

03 

HP 

a 

o 

CO 

'hp 

o 

4 

co 

bO 

.2 

*d 

d 

H 

*n 

HP 

CO 

O 

CO 

’ CO 

d 

0 

u 

• rH 

"d 

a 

simulations  for  October  6 - 11,  1987.  The  models  were  run  at  a 0.5  hour  time  step  to 
capture  the  diurnal  cycle  in  potential  evapotranspiration.  Three  models  are  compared:  a 
fully  distributed  model  in  which  the  spatial  patterns  of  the  soil-topographic  parameters 
and  rainfall  are  preserved;  a macroscale,  distributed  model  in  which  the  spatial  variability 
in  the  soil-topographic  parameters  is  accounted  for  statistically;  and  an  aggregated,  one- 


dimensional  model  in  which  parameters  and  inputs  are  spatially  constant  at  their  average 


In  the  previous  section  through  a discussion  of  Figure  1,  a qualitative  analysis  was  pre- 
sented concerning  the  mechanisms  leading  to  a non-linear  behavior  of  evaporation  and  the 
under  estimation  of  the  evaporative  flux  when  aggregation  is  carried  out.  Close  examina- 
tion of  the  concepts  around  Figure  1 leads  to  the  insight  that  difficulties  in  aggregation 


simulation,  yet  still  preserve  the  essential  features  of  the  aggregation  problem. 


soil  moisture  less  than  9 * is  defined  as  F,  and  will  be  referred  to  as  the  traction  ot  tne 
watershed  under  ’soil-controlled’  evaporation  since  the  evaporation  from  these  areas  is 
limited  by  the  soil- vegetation  system.  The  fraction  of  the  watershed  1 — F,  is  under 
’atmospheric-controlled’  evaporation  and  the  evaporation  rate  is  at  the  potential  rate. 

Let  the  ratio  of  actual  evaporation  to  potential  evaporation  for  soil  moisture,  9„  less 


than  9*  be  represented  as 
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3.2  Derivation  of  the  Evaporative  Fraction. 
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The  sensitivity  of  the  evaporation  scaling  ratio  depends  on  three  parameters:  the  ratio  ot 
the  equilibrium  soil  moisture  0*to  the  mean  soil  moisture  (9),  8*/ (8)]  the  soil- vegetation 


sensitivity  factor  to  soil  moisture  stress,  /3,and  the  variability  of  soil  moisture  within  the 
region  measured  by  its  coefficient  of  variation,  Cv. 

As  discussed  earlier,  9 * represents  the  soil  moisture  level  at  which  the  soil-vegetation 
system  is  in  equilibrium  with  the  atmospheric  evaporative  demand.  Thus,  areas  with  soil 
moisture  lower  than  8*  will  exhibit  a soil- vegetation  control  and  have  an  actual  evaporative 
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Shuttle  Imaging  Radar  soil  moisture  experiment  in  the  same  basin,  field  (0.25  to  0.5  sq. 
mi.)  observed  soil  moisture  Cv  (from  approximately  20  to  25  samples  per  field)  ranged 
from  0.11  (for  winter  wheat  in  dry  conditions)  to  0.44  for  natural  pasture  after  a light 
rain.  The  five  winter  wheat  fields  ranged  from  0.11  to  0.24,  the  five  pasture  fields  from 
0.15  to  0.44  and  two  bare  soil  fields  (plowed)  from  0.19  to  0.28.  The  fields  have  uniform 


vegetation  type  and.  soil  texture  class. 

A series  of  sensitivity  experiments  were  conducted  by  varying  these  three  parameters: 
9*/ (9)  within  the  range  of  0.2  to  2;  /3  within  the  range  0.05  to  0.6;  and  Cv  Within  the 
range  0.10  to  0.40.  The  results  are  provided  in  Figures  3 and  4.  Figure  3 provides  the 
sensitivity  to  changes  in  /3  for  different  Cv  values,  while  Figure  4 provides  the  sensitivity 
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during  mid-day  Summer  periods  when  atmospheric  evaporative  demands  are  high  and 
surface  soil  moisture  low.  Results  indicate  that  the  aggregated  model  underestimates 
evaporation  in  the  range  of  15%  to  50%  relative  to  the  modeled  evaporation  from  the 
distributed  model  for  similar  conditions. 

(iii)  For  cases  of  low  (3  and  low  variation  in  soil  moisture,  the  evaporative  scaling 


coefficient  is  close  to  1.0,  showing  that  the  aggregated  models  and  the  distributed  models 
estimate  the  same  evaporation. 

These  results  imply  that  aggregated,  lumped  parameter  models  will  rarely  compute 
evaporation  that  is  equal  to  that  computed  by  a distributed  model;  even  under  appar- 
ently ’atmospheric  controlled’  conditions.  This  conclusion  is  due  to  the  variability  m soil 
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Dubayah  et  al.,  (1996)  performed  a statistical  self-similarity  analysis  using  soil  mois- 
ture fields  estimated  from  an  airborne  remote  sensing  passive  microwave  (ESTAR)  in- 

s 

strument.  The  data  was  collected  over  the  Little  Washita,  OK  watershed  during  June 
1992.  Results  show  that  the  soil  moisture  fields  scale,  but  in  a non-linear  manner;  i.e. 
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Figure  la:  Schematic  of  deriving  the  probability  distribution  for  evaporation  from  the 
distribution  of  water  table  depths  for  dry  conditions  (average  water  table  depth  is  3.5 
m).  (i)  distribution  of  the  water  table  depths,  (ii)  surface  soil  moisture  as  a function  of 
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